Transverse coherence properties of X-ray beams in 
third-generation synchrotron radiation sources 



Gianluca Geloni Evgeni Saldin Evgeni Schneidmiller 
and Mikhail Yurkov 

Deutsches Elektronen-Synchrotron (DESY), Hamburg, Germany 



Abstract 

This article describes a complete theory of spatial coherence for undulator ra- 
diation sources. Current estimations of coherence properties often assume that 
undulator sources are quasi-homogeneous, like thermal sources, and rely on the 
application of the van Cittert-Zernike theorem for calculating the degree of trans- 
verse coherence. Such assumption is not adequate when treating third generation 
light sources, because the vertical (geometrical) emittance of the electron beam is 
comparable or even much smaller than the radiation wavelength in a very wide 
spectral interval that spans over four orders of magnitude (from 0.1 A up to lO'^A). 
Sometimes, the so-called Gaussian-Schell model, that is widely used in statistical 
optics in the description of partially-coherent sources, is applied as an alternative to 
the quasi-homogeneous model. However, as we will demonstrate, this model fails 
to properly describe coherent properties of X-ray beams from non-homogeneous 
undulator sources. As a result, a more rigorous analysis is required. We propose 
a technique, based on statistical optics and Fourier optics, to explicitly calculate 
the cross-spectral density of an undulator source in the most general case, at any 
position after the undulator. Our theory, that makes consistent use of dimension- 
less analysis, allows relatively easy treatment and physical understanding of many 
asymptotes of the parameter space, together with their region of applicability. Par- 
ticular emphasis is given to the asymptotic situation when the horizontal emittance 
is much larger than the radiation wavelength, and the vertical emittance is arbi- 
trary. This case is practically relevant for third generation synchrotron radiation 
sources. 
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1 Introduction 



In recent years, continuous evolution of synchrotron radiation (SR) sources 
has resulted in a dramatic increase of brilliance with respect to older designs. 
Among the most exciting properties of third generation facilities of today is 
a high flux of coherent X-rays [1]. The availability of intense coherent X-ray 
beams has triggered the development of a number of new experimental 
techniques based on coherence properties of light such as scattering of 
coherent X-ray radiation. X-ray photon correlation spectroscopy and phase 
contrast imaging. The interested reader may find an extensive reference list 
to these applications in [2]. 

Characterization of transverse coherence properties of SR at a specimen 
position is fundamental for properly planning, conducting and analyzing 
experiments involving above-mentioned techniques. These objectives can 
only be met once coherence properties of undulator sources are character- 
ized and then propagated along the photon beamline up to the specimen. 

This paper is mainly dedicated to a description of transverse coherence 
properties of SR sources. Therefore, it constitutes the first step for tracking 
coherence properties through optical elements, which can be only done 
when the source is quantitatively described. Solution of the problem of 
evolution of radiation properties in free-space is also given. This can be 
used to characterize radiation at the specimen position when there are no 
optical elements between the source and the specimen. 

Since SR is a random process, the description of transverse coherence proper- 
ties of the source and its evolution should be treated in terms of probabilistic 
statements. Statistical optics presents most convenient tools to deal 

with fluctuating electromagnetic fields. However, it was mainly developed 
in connection with polarized thermal light that is characterized by quite 
specific properties. Besides obeying Gaussian statistics, polarized thermal 
light has two other specific characteristics allowing for major simplifications 
of the theory: stationarity and homogeneity (i.e. perfect incoherence) of the 
source 

As we will see, SR obeys Gaussian statistics, but stationarity and homo- 
geneity do not belong to SR fields. Thus, although the language of statistical 
optics must be used to describe SR sources, one should avoid a-priori in- 
troduction of an incorrect model. In contrast to this, up to now it has been 
a widespread practice to assume that undulator sources are perfectly inco- 
herent (homogeneous) and to draw conclusions about transverse coherence 
properties of undulator light based on these assumptions [|2fl. 

In particular, in the case of thermal light, transverse coherence properties 
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of the radiation can be found with the help of the well-known van Citter- 
Zernike (VCZ) theorem. However, for third generation light sources either 
planned or in operation, the horizontal electron beam (geometrical) emit- 
tanc^ is of order of 1 -r 3 nm while the vertical emittance is bound to the 
horizontal through a coupling factor C ~ 10"^, corresponding to vertical 
emittances of order 0.1 -r 0.3 A. As we will show, these facts imply that 
the VCZ theorem cannot be applied in the vertical direction, aside for the 
hard X-ray limit at wavelengths shorter than 0.1 A. Similar remarks hold 
for future sources, like Energy Recovery Ltnac (ERL)-based spontaneous 
radiators. ERL technology is expected to constitute the natural evolution 
of today SR sources and to provide nearly fully diffraction-limited sources 
in the lA-range, capable of three order of magnitudes larger coherent flux 
compared to third generation light sources [1]. Horizontal and vertical emit- 
tance will be of order of 0.3A, which rules out the applicability of the VCZ 
theorem. 

The need for characterization of partially coherent undulator sources is em- 
phasized very clearly in reference [6]. However, studying spatial-coherence 
properties of undulator sources is not a trivial task. Difficulties arise when 
one tries to simultaneously include the effect of intrinsic divergence of un- 
dulator radiation, and of electron beam size and divergence. 

An attempt to find the region of applicability of the VCZ theorem for third 
generation light sources is reported in [7]. Authors of |0] conclude that 
the VCZ theorem can only be applied when SR sources are close to the 
diffraction limit. We will see that the VCZ theorem can only be applied in 
the opposite case. 

A model to describe partially-coherent sources, the Gaussian-Schell model, 
is widely used in statistical optics. Its application to SR is described in 
iB&lloll- The Gaussian-Schell model includes non-homogeneous sources, 
because the typical transverse dimension of the source can be comparable or 
smaller than the transverse coherence length. While are of general 

theoretical interest, they do not provide a satisfactory approximation of 
third generation SR sources. In fact, as we will see, undulator radiation has 
very specific properties that cannot be described in terms of Gaussian-Schell 
model. 

The first treatment of transverse coherence properties from SR sources ac- 
counting for the specific nature of undulator radiation and anticipating 
operation of third generation SR facilities is given, in terms of Wigner dis- 



Emittances are measured in m • rad. However, radians are a dimensionless unit. 
Therefore, we can present emittances measured in meters. This is particularly useful 
in the present paper, since we need to compare emittances with radiation wave- 
length. 
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tribution, in l^]. References lU, 12] present the most general algorithm 
for calculating properties of undulator radiation sources. 



In the present paper we base our theory on the characterization of the 
cross-spectral density of the system. The cross-spectral density is merely 
the statistical correlation function of the radiation field at two different po- 
sitions on the observation plane at a given frequency. It is equivalent to the 
Wigner distribution, the two quantities being related by Fourier transforma- 
tion. Based on cross-spectral density we developed a comprehensive theory 
of third generation SR sources in the space-frequency domain, where we 
exploited the presence of small or large parameters intrinsic in the descrip- 
tion of the system. First, we took advantage of the particular but important 
situation of perfect resonance, when the field from the undulator source 
can be presented in terms of analytical functions. Second, we exploited 
the practical case of a Gaussian electron beam, allowing further analytical 
simplifications. Third, we took advantage of the large horizontal emittance 
(compared with the radiation wavelength) of the electron beam, allowing 
separate treatments of horizontal and vertical directions. Finally, we stud- 
ied asymptotic cases of our theory for third generation light sources with 
their region of applicability. In particular, we considered both asymptotes 
of a small and a large vertical emittance compared with the radiation wave- 
length, finding surprising results. In the case of a small vertical emittance 
(smaller than the radiation wavelength) we found that the radiation is not 
diffraction limited in the vertical direction, an effect that can be ascribed 
to the influence of the large horizontal emittance on the vertical coherence 
properties of radiation. In the case of a large vertical emittance we described 
the quasi-homogeneous source in terms of a non-Gaussian model, which has 
higher accuracy and wider region of applicability with respect to a simpler 
Gaussian model. In fact, the non-Gaussian model proposed in our work ac- 
counts for diffraction effects. Nonetheless, in order to solve the propagation 
problem in free-space, a geometrical optics approach can still be used. 



We organize our work as follows. First, a general algorithm for the calcu- 
lation of the field correlation function of an undulator source is given in 
Section 2. Such algorithm is not limited to the description of third genera- 
tion light sources, but can also be applied, for example, to ERL sources. In 
Section 3 we specialize our treatment to third generation light sources ex- 
ploiting a large horizontal emittance compared to the radiation wavelength, 
and studying the asymptotic limit for a small vertical emittance of the elec- 
tron beam. In the following Section 4 we consider the particular case of 
quasi-homogeneous undulator sources, and discuss the applicability range 
of the van VCZ theorem. Finally, in Section 5 we come to conclusions. 



4 



thermal 
light 



completely chaotic 
polarized light 



lamp 



polarizer 



spatial 
filter 



spectral 
filter 



[(r,t) 



stationary Gaussian 
random process 



time-averaging 
intensity 




t 

Fig. 1. Light intensity from an incandescent lamp driven by a constant electric 
current. A statistically stationary wave has an average that does not vary with 
time. 

2 Cross-spectral density of an undulator source 



2.1 Thermal light and Synchrotron Radiation: some concepts and definitions 



The great majority of optical sources emits thermal light. Such is the case of 
the sun and the other stars, as well as of incandescent lamps. This kind of 
radiation consists of a large number of independent contributions (radiating 
atoms) and is characterized by random amplitudes and phases in space and 
time. Its electromagnetic field can be conveniently described in terms of 
statistical optics that has been intensively developed during the last few 
decades (siH. 

Consider the light emitted by a thermal source passing through a polariza- 
tion analyzer, as in Fig. 1. Properties of polarized thermal light are well- 
known in statistical optics, and are referred to as properties of completely 
chaotic, polarized light 10, H]- Thermal light is a statistical random process 
and statements about such process are probabilistic statements. Statistical 
processes are handled using the concept of statistical ensemble, drawn from 
statistical mechanics, and statistical averages are performed over an ensem- 
ble, or many realizations, or outcomes of the statistical process under study. 
Polarized thermal sources obey a very particular kind of random process in 
that it is Gaussian and stationar}|il Moreover, they are homogeneous. 

^ Ergodic too. From a qualitative viewpoint, a given random process is ergodic 
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The properties of Gaussian random processes are well-known in statistical 
optics. For instance, the real and imaginary part of the complex amplitudes 
of the electric field from a polarized thermal source have Gaussian distribu- 
tion, while the instantaneous radiation power fluctuates in accordance with 
the negative exponential distribution. It can be shown ^ that a linearly 
filtered Gaussian process is also a Gaussian random process. As a result, the 
presence of a monochromator and a spatial filter as in the system depicted 
in Fig. 1 do not change the statistics of the signal. Finally, higher order field 
correlation functions can be found in terms of the second order correlation. 
This dramatically simplifies the description of the random process. 

Stationarity is a subtle concept. There are different kinds of stationarity. 
However, for Gaussian processes different kinds of stationarity coincide 
10, 0]- In this case, stationarity means that all ensemble averages are inde- 
pendent of time. It follows that a necessary condition for a certain process to 
be stationary is that the signal last forever. Yet, if a signal lasts much longer 
than the short-scale duration of the field fluctuations (its coherence time t^) 
and it is observed for a time much shorter than its duration Oj, but much 
longer than Tc it can be reasonably considered as everlasting and it has a 
chance to be stationary as well, as in the case of thermal light. 

Finally, thermal sources are homogeneous. The field is correlated on the 
transverse direction over the possible minimal distance, which is of order of 
a wavelength. This means that the radiation intensity at the source remains 
practically unvaried on the scale of a correlation length. One can say, equiv- 
alently, that the source is homogeneous. Homogeneity is the equivalent, in 
the transverse direction, of stationarity and implies a constant ensemble- 
averaged intensity along the transverse direction. 

Consider now a SR source, as depictec^Elin Fig. 2. Like thermal light, also SR 
is a random process. In fact, relativistic electrons in a storage ring emit SR 
passing through bending magnets or undulators. The electron beam shot 
noise causes fluctuations of the beam density which are random in time and 
space from bunch to bunch. As a result, the radiation produced has random 
amplitudes and phases. Moreover, in Section 2.2 we will demonstrate that 
SR fields obey Gaussian statistics. Statistical properties satisfied by single- 



when all ensemble averages can be substituted by time averages. Ergodicity is a 
stronger requirement than stationarity ii. 

^ Radiation at the detector consists of a carrier modulation of frequency co subjected 
to random amplitude and phase modulation. The Fourier decomposition of the 
radiation contains frequencies spread about the monochromator bandwidth: it is 
not possible, in practice, to resolve the oscillations of the radiation fields which 
occur at the frequency of the carrier modulation. Therefore, for comparison with 
experimental results, we average the theoretical results over a cycle of oscillations 
of the carrier modulation. 
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Fig. 2. The intensity of an X-ray beam from a SR source. A statistically non-stationary 
wave has a time-varying intensity averaged over an ensemble of bunches. 



electron contributions (elementary phasors) to the total SR field are weaker 
than those satisfied by single-atom contributions in thermal sources. Thus, 
the demonstration that thermal light obeys Gaussian statistics cannot be 
directly applied to the SR case, and some condition should be formulated 
on the parameter space to define the region where SR is indeed a Gaussian 
random process. We will show this fact with the help of Appendix A of [IS,] . 

In contrast with thermal light, SR is intrinsically non-stationary, because it 
presents a time-varying ensemble-averaged intensity on the temporal scale 
of the duration of the X-ray pulse generated by a single electron bunch. For 
this reason, in what follows the averaging brackets (...) will always indicate 
the ensemble average over bunches (and not a time average). Finally, SR 
sources are not completely incoherent, or homogeneous. In fact, there is a 
close connection between the state of coherence of the source and the angular 
distribution of the radiant intensity. A thermal source that is correlated 
over the minimal possible distance (which is of order of the wavelength) is 
characterized by a radiant intensity distributed over a solid angle of order 
271. This is not the case of SR light, that is confined within a narrow cone in 
the forward direction. The high directionality of SR rules out the possibility 
of description in terms of thermal light. As we will see, depending on the 
situation, SR may or may not be described by a quasi-homogeneous model, 
where sources are only locally coherent over a distance of many wavelengths 
but the linear dimension of the source is much larger than the correlation 
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distance0- 



In spite of differences with respect to the simpler case of thermal light, SR 
fields can be described in terms of statistical optics. However, statistical 
optics was developed in relation with thermal light emission. Major as- 
sumptions typical of this kind of radiation like, for example, stationarity, 
are retained in textbooks [3, 4]. Thus, usual statistical optics treatment must 
be modified in order to deal with SR problems. We will reduce the problem 
of characterization of transverse coherence properties of undulator sources 
in the space-frequency domain to the calculation of the correlation of the 
field produced by a single electron with itself. This correlation is known 
as cross-spectral density. The non-stationarity of the process imposes some 
(practically non restrictive) condition on the parameter space region where 
a treatment based on cross-spectral density can be applied to SR phenom- 
ena. Such condition involves the length of the electron bunch, the number 
of undulator periods and the radiation wavelength. Given the fact that the 
electron bunch length can vary from about 10 ps for third generation SR 
sources to about 100 fs for ERL sources, understanding of the region of ap- 
plicability this condition is a-priori not obvious. All these subjects will be 
treated in the next Section 2.2. 



2.2 Second-order correlations in space-frequency domain 

In SR experiments with third generation light sources, detectors are limited 
to about 100 ps time resolution. Therefore, they cannot resolve a single 
X-ray pulse in time domain, whose duration is about 30 ps. They work, 
instead, by counting the number of photons at a certain frequency over 
an integration time longer than the pulse. It is therefore quite natural to 
consider signals in the frequency domain. With this in mind, let Ei,{z,r,co) 
be a fixed polarization component of the Fourier transform of the electric 
field at location (z, r), in some cartesian coordinate system, and frequency 
cu by a given collection of electromagnetic sources. We will often name 
it, slightly improperly, "the field". Subscript "h" indicates that the field is 
generated by the entire bunch. Ei,{z,r,co) is linked to the time domain field 
£b(z, r, t) through the Fourier transform 



Eh{co) = I dtEi,{t) exp{icot) , Ei,{t) = — I dcoEi,{co) exp{-icot) . (1) 



^ Note that the high directionality of SR is not in contrast with the poor coherence 
which characterizes the quasi-homogeneous limit. 
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We will be interested in the case of an ultra relativistic electron beam going 
through a certain magnetic system, an undulator in particular. In this case z 
is the observation distance along the optical axis of the undulator and r fixes 
the transverse position of the observer. The contribution of the k-th electron 

to the field depends on the transverse offset k and deflection angles rjk that 
the electron has at some reference point on the optical axis z, e.g. the center 
of the undulator. Moreover, the arrival time tk at the center of the undulator 
has the effect of multiplying the field by a phase factor exp {iivtk), i.e. the 
time-domain electric field is retarded by a time tk. At this point we do not 
need to explicitly specify the dependence on offset and deflection. The total 
field can be written as 

Eb (z, r,a)) = ^E (^ifk, Ik, z, r, co) exp (icotk) , (2) 

/c=l 



where ifk, h and are random variables and Ne is the number of electrons in 
the bunch. Note that E in Eq. (2) is a complex quantity, that can be written 
as E -\ Ek \ exp{i(pk)- It follows that the SR field pulse at fixed frequency 
and position is a sum of a many phasors, one for each electron, of the form 
E exp (icotk) =1 Ek \ exp{i(pk + icotk). 

Elementary phasors composing the sum obey three statistical properties, 
that are satisfied in SR problems of interest. First, random variables tk are 

statistically independent of each other, and of variables rjk and k- Second, 
random variables | Ek \ (at fixed frequency co), are identically distributed for 
all values of k, with a finite mean {| Ek |) and a finite second moment {| Ek p). 
These two assumptions follows from the properties of shot noise in a storage 
ring, which is a fundamental effect related with quantum fluctuations. Third, 
we assume that the electron bunch duration oj is large enough so that 
coot » 1: under this assumption, phases cotk can be regarded as uniformly 
distributed on the interval {0,2n). The assumption cooj » 1 exploits the 
first large parameter of our theory, and is justified by the fact that co is the 
undulator resonant frequency, which is high enough, in practical cases of 
interest, to guarantee that coaj » 1 for any realistic choice of aj- Based on 
the three previously discussed properties, and with the help of the central 
limit theorem, it can be demonstratec^l] that the real and the imaginary part 
of Efc are distributed in accordance to a Gaussian law. 



^ The proof follows from a slight generalization of Section 2.9 in Namely, it can 
be shown by direct calculation that real and imaginary part of the total phasor 
are uncorrelated, with zero mean and equal variance. Then, using the central limit 
theorem, we conclude that is a circular complex Gaussian random variable (at 
fixed z, r and co). 
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As a result, SR is a non-stationary Gaussian random process. 

Because of this, higher-order correlation functions can be expressed in terms 
of second-order correlation functions with the help of the moment theorem 
1^. As a result, the knowledge of the second-order field correlation function 
in frequency domain, T^,, is all we need to completely characterize the signal 
from a statistical viewpoint. The following definition holds: 

To; (z, ri, r2, cvi, 0)2) = {Eh (z, ri, cvi) El (z, r2, (V2)) , (3) 

where brackets < ... > indicate ensemble average over electron bunches. For 
any given function w (rfk, h, tk), the ensemble average is defined as 



00 

{iv(r]k,h,tk)) = J drfk J dk J dtkw{fikJk,tk)p(r]k,h,tk) , (4) 



where integrals in dh and drjk span over all offsets and deflections, and P = 
P{T]k, h, h) indicates the probability density distribution in the joint random 
variables -qk, h, and h. Note that, since all electrons have the same probability 
of arrival around a given offset, deflection, and time, P is independent of 

k. Moreover, already discussed independence of i)c from 1;^ and allows to 
write P as 



P Ik, tk) = fi. (5, fik) f{tk) . (5) 

Here / is the longitudinal bunch profile of the electron beam, while fj_ is the 
transverse phase-space distribution. 

Substituting Eq. (2) in Eq. (3) one has 



r^, = ( 2], E(']m,Zm,z,j^i,ti^i)E*(?]n,Z„,z,r2,a;2)exp[zXa;i^„, - a;2^„)])- 

(6) 

Expansion of Eq. (6) gives 



(e (rf„„ z, ri, cui) E* {i]„„ /„„ z, ^2, exp [i{a)i - C02)t„,]) 

m=l 
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+ 2], (?fm, Im, z. A, 00^ exp (za;ii„,)) (e* (if„, /„, z, r2, 0^2) exp {-ia)2t„)) . 



(7) 



With the help of Eq. (4) and Eq. (5), the ensemble average (exp (icotk)) can be 
written as the Fourier transform of the longitudinal bunch profile function 
/, that is 



00 

{exp (icotk)) = J dtkf{tk)exp{icotk) = f{co) 



Using Eq. (8), Eq. (7) can be written as 



(8) 



= /(cui - 0)2) (e (rf,„, z, Ti, cui) E* (rf„, /„, z, r2, ct^i)) 

m=l 

+ f{(0i)f{-(02) (e z, ri, a;i)) (e* (?]„, C ^2, 0^2)) , (9) 



where f*{oJ2) = f{-0J2) because / is a real function. When the radiation 
wavelengths of interest are much shorter than the bunch length we can 
safely neglect the second term on the right hand side of Eq. (9) since the 
form factor product f{a)i)f {-0)2} goes rapidly to zero for frequencies larger 
than the characteristic frequency associated with the bunch length: think 
for instance, at a centimeter long bunch compared with radiation in the 
Angstrom wavelength rang^. Therefore we write 



Ta, = ^ [{cOi - CO2) (E (rfm, Im, Z, A, COlj E* (ifm, Im, Z, A, (^l)) 

m=l 

= Ne f{a)i - 0)2) (e (rf, I, z, A, f^i) E* (rf, I, z, A, oJi^ ■ (10) 

As one can see from Eq. (10) each electron is correlated just with itself: cross- 
correlation terms between different electrons was included in the second 
term on the right hand side of Eq. (9), that has been dropped. 

^ When the radiation wavelength of interested is comparable or longer than the 
bunch length, the second term in Eq. (9) is dominant with respect to the first, 
because it scales with the number of particles squared: in this case, analysis of the 
second term leads to a treatment of coherent synchrotron radiation phenomena 
(CSR). In this paper we will not be concerned with CSR and we will neglect the 
second term in Eq. (9), assuming that the radiation wavelength of interest is shorter 
than the bunch length. Also note that /(cui - (02) depends on the difference between 
(x)\ and CU2, and the first term cannot be neglected. 



11 



- 10" -10 



-3 



G(i^, 






f(cj -co') 




^^10" for hard X-rays 



Fig. 3. Schematic illustration of the relative frequency dependence of the spec- 
tral correlation function f{a) — co') and of the cross-spectral density function (the 
cross-power spectrum) G(z, fi, A, (o) of SR at points ri and A at frequency co. 

If the dependence of E on cu is slow enough, E does not vary appreciably 
on the characteristic scale of /. Thus, we can substitute E*{rj, I, z, r2, CO2) with 
E*{ri,l,z,r2,coi) in Eq. (10). The situation is depicted in Fig. 3. On the one 
hand, the characteristic scale of / is given by 1 / Oj, where oj is the character- 
istic bunch duration. On the other hand, the minimal possible bandwidth of 
undulator radiation is achieved on axis and in the case of a filament beam. It 
is peaked around the resonant frequency co,- = 2y^c/A-^ (A-^ being the undu- 
lator period and the average longitudinal Lorentz factor) and amounts to 
<^r/Nw> ^ 1 being the number of undulator periods (of order 10^ - lO''). 
Since C0r/Nj„ is a minimum for the radiation bandwidth, it should be com- 
pared with 1 /ctt. For instance, at wavelengths of order lA, N^, ~ 10^ and 
30 ps (see ml 

one has cOr/Nw ~ 2 • lO^'' Hz, which is much larger than 
l/cTj ~ 3 • 10^" Hz. From this discussion follows that, in practical situations 
of interest, we can simplify Eq. (10) to 

r^,(z, ri, r2, coi, C02) = N, f{coi - C02)G (z, A, coi) , (11) 



Eq. (11) fully characterizes the system under study from a statistical view- 
point. Correlation in frequency and space are expressed by two separate 
factors. In particular, spatial correlation is expressed by the cross-spectral 



where 




(12) 
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Fig. 4. Measurement of the cross-spectral density of a undulator source, (a) Young's 
double-pinhole interferometer demonstrating the coherence properties of undula- 
tor radiation. The radiation beyond the pinholes must be spectrally filtered by a 
monochromator or detector (not shown in figure), (b) The fringe visibility of the 
resultant interference pattern. 

density functiorH] G. In other words, we are able to deal separately with spa- 
tial and spectral part of the correlation function in space-frequency domain 
under the non-restrictive assumption (OrlNw » lloj. 

Eq. (1 1) is the result of our theoretical analysis of the second-order correlation 
function in the space-frequency domain. We can readily extend this analysis 
to the case when the observation plane is behind a monochromator with 
transfer function T{a)). In this case, Eq. (11) modifies to 



(z, ri, r2, coi, (1)2) = Neficoi - C02)G (z, n, r2, coi) T{coi)T{co2) . 



(13) 



It is worth to note that, similarly to Eq. (11), also Eq. (13) exhibits separability 
of correlation functions in frequency and space. 

From now on we will be concerned with the calculation of the cross-spectral 
density G(z, ti, Y2, oS), independently of the shape of the remaining factors 
on the right hand side of Eq. (13). 



Note, however, that G depends on co. 
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Before proceeding we introduce, for future reference, the notion of spectral 
degree of coherence g, that can be presented as a function of ri and rz: 

(|£«)f) (|E(r-,)f^ 



Consider Fig. 4, depicting a Young's double-pinhole interferometric mea- 
sure. Results of Young's experiment vary with and r2- The modulus of the 
spectral degree of coherence, \g\ is related with the fringe visibility of the 
interference pattern. In particular, the relation between fringes visibility V 
and g is given by 



/ 1/2 / t\ 1/2 

E(n)| )+ |E(rl)| 



Phase of g is related to the position of the fringes instead. Thus, spectral 
degree of coherence and cross-spectral density are related with both ampli- 
tude and position of the fringes and are physically measurable quantities 
that can be recovered from a Young's interference experiment. 



2.3 Relation between space-frequency and space-time domain 

This paper deals with transverse coherence properties of SR sources in 
the space-frequency domain. However, it is interesting to briefly discuss 
relations with the space-time domain, and concepts like quasi-stationarity, 
cross-spectralpurity and quasi-monochromaticity that are often considered 
in literature iSi]. 

First, the knowledge of F^ in frequency domain is completely equivalent to 
the knowledge of the mutual coherence function fl^ : 

rt(z. A, ri, ti, tz) = {E{z, A, ti)E*{z, A, ii)) ■ (16) 
Next to Ft, the complex degree of coherence is defined as 



yt{z, ri, r2, ti, t2) = — — -j^ . (17) 

[Ff (z, ri, ri, ti, ti)Tt{z, r2, ri, tz, tz)] 
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The presence of a monochromator (see Eq. (13)) is related with a bandwidth 
of interest AcUm, centered around a given frequency cOg (typically, the undu- 
lator resonant frequency). Then, T{ijj) is peaked around cuo and rapidly goes 
to zero as we move out of the range {cOo - Aco^/2,cOo + AcOm/^). Using Eq. 
(13) we write the mutual coherence function as 



^t{z,n,r2,h,h)^-^^ J dcoi J dco2f{coi - co2)T{coi)r{co2) 

—oo —oo 

xG(z, ri, r2, (02) exp {-icoiti) exp {ico2t2) • (18) 

If the characteristic bandwidth of the monochromator, Acora, is large enough 
so that T does not vary appreciably on the characteristic scale of /, i.e. 
Acum » I/ct/ then T{coi)T*{co2)f{coi - 0)2) is peaked at cui - 0)2 = 0. In 
this case the process is quasi-stationary With the help of new variables 
Aoj -0)1-0)2 and o) = (oji + 0)2) /2, we can simplify Eq. (18) accounting 
for the fact that f{a)i - 0)2) is strongly peaked around Ao) = 0. In fact we 
can consider T{o)i)T*{o)2)G{z,ri,r2,<^i) - |T((D)pG(z,ri,r2/(^>)/ so that we can 
integrate separately in Acu and d) to obtain 

00 

Tf (z, ri, f2, h, t2) = -^A^ J dAo) f{Ao)) exp {-iAo)t) 

—00 

00 

X Jdo) |T((D)pG(z, ri, ^2, d)) exp [-icoAt] . (19) 

— CX) 

where t = (ti + 12)/ 2 and At = h - t2. This means that the mutual co- 
herence function can be factorized as rt(F, At) = f{t)Gt{At). In the case no 
monochromator is present, Gf coincides with the Fourier transform of the 
cross-spectral density, and the correspondent correlation function Tc^{a}, Ao)) 
has been seen to obey Eq. (10). Therefore: 

(fli) The temporal correlation function Gt{At) and the spectral density distri- 
bution of the source \T{o))\^G{cD) form a Fourier pair. 

(&i) The intensity distribution of the radiation pulse f{t) and the spectral 
correlation function / {Aa>) form a Fourier pair. 

Statement (fli) can be regarded as an analogue, for quasi-stationary sources, 
of the well-known Wiener-Khinchin theorem, which applies to stationary 
sources and states that the temporal correlation function and the spectral 
density are a Fourier pair. Since there is symmetry between time and fre- 
quency domains, an inverse Wiener-Khinchin theorem must also hold, and 
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Fig. 5. Reciprocal width relations of Fourier trarisform pairs. 

can be obtained by the usual Wiener-Khinchin theorem by exchanging fre- 
quencies and times. This is simply constituted by statement (fci). Intuitively, 
statements (fli) and (bi) have their justification in the reciprocal width rela- 
tions of Fourier transform pairs (see Fig. 5). 

It should be stressed that statistical optics is almost always applied in the 
stationary case f^,^. Definitions of Ff and yt are also usually restricted to 
such case. The case of a stationary process corresponds to the asymptote of a 
Dirac 6-function for the spectral correlation function /. The inverse Wiener- 
Khinchin theorem applied to this asymptotic case would imply an infinitely 
long radiation pulse, i.e. an infinitely long electron beam. In contrast to this, 
the width of / is finite, and corresponds to a finite width of / of about 30 
ps. Thus, one cannot talk about stationarity. However, when AcOm ^ I/ct/ 
the spectral width of the process (i.e. the width of \T\^G in co) is much larger 
than the width of /. In this case, the process is quasi-stationary. The situation 
changes completely if a monochromator with a bandwidth AcUm ^ 1 /oj is 
present. In this case, Eq. (19) cannot be used anymore, and one is not allowed 
to treat the process as quasi-stationary. In the large majority of the cases 
monochromator characteristics are not good enough to allow resolution of 
/. There are, however, some special cases when AcOm ^ I/ct- For instance. 



in llSy a particular monochromator is described with a relative resolution 
of 10"^ at wavelengths of about lA, or a;,, ~ 2 • 10^^ Hz. Let us consider, as 
in lE^, the case of radiation pulses of 32 ps duration. Under the already 
accepted assumption l/aj ^ cUo/Na,, we can identify the radiation pulse 
duration with oj- Then we have Aco^ ~ 2 • 10^^ Hz which is of order of 
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InjoT ~ 2 • 10^^ Hz: this means that the monochromator has the capability 
of resolving /. 

Cases discussed up to now deal with radiation that is not cross-spectrally 
pure |[l6l] (or paragraph 4.5.1). In fact, the absolute value of the spectral 
degree of coherence g is a function of co. Moreover, as remarked in ||8[l9l,[lol], 
the spectrum of undulator radiation depends on the observation point. This 
fact can also be seen in the time domain from Eq. (18) or Eq. (19), because the 
complex degree of coherence yt cannot be split into a product of temporal 
and spatial factors. However, if we assume AojmNiv/cOo ^ 1 (that is usually 
true), G(z, ri, r2, o)) = G{z, Ti, Ti, oJq) is a constant function of frequency within 
the monochromator line, i.e. it is independent on the frequency co. As a 
result Ff in Eq. (18) can be split in the product of a temporal and a spatial 
factor and therefore, in this case, light is cross-spectrally pure: 

Tf (Zo, ToX, Tol, h, tl) = Nejtih, ^2)G (z, fi, fi, COo) , (20) 



where gt{ti, ti) is defined by 



oo oo 

St = j^^ J ^^^1 J dco2f{coi-co2)T{coi)r{co2)exp[i{co2t2-coiti)] 



(21) 



Note that, for instance, in the example considered before Acom/coo - 10 ^ 
and N^u - 10^ 10^, i.e. AcOj^^Nra/cOo « 1. 

It is important to remark that, since we are dealing with the process in the 
space-frequency domain, whether the light is cross-spectrally pure or not 
is irrelevant concerning the applicability of our treatment, because we can 
study the cross-spectral density G for any frequency component. 

Finally, for the sake of completeness, it is interesting to discuss the relation 
between G and the mutual intensity function as usually defined in text- 
books Si in quasimonochromatic conditions. The assumption Acum » l/crj 
describes a quasi-stationary process. In the limit Ot — > we have a sta- 
tionary process. Now letting AcOm — > slowly enough so that Ao^m » I/ct/ 
Eq. (19) remains valid while both /(Acu) and |T((D)p become approximated 
better and better by Dirac 6-functions, b{Aa)) and 6(d) - cOo), respectively. 
Then, Ft ~ G exp[-za;o(ii - ^2)]- Aside for an unessential factor, depending 
on the normalization of / and |T(a;)p, this relation between Ff and G allows 
identification of the mutual intensity function with G as in [S^^. In this case, 
light is obviously cross-spectrally pure. 
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2.4 Undulator field by a single particle with offset and deflection 

In order to give an explicit expression for the cross-spectral density of un- 
dulator radiation, we first need an explicit expression for E [f], I, z, r, a;), the 
field contribution from a single electron with given offset and deflection. 
This can be obtained by solving paraxial Maxwell's equations in the space- 
frequency domain with a Green's function technique. We refer the reader 
to [17J, where it was also shown that paraxial approximation always ap- 
plies for ultrarelativistic systems with « 1, y being the relativistic 
Lorentz factor. Since paraxial approximation applies, the envelope of the 
field £((!;) = £(a;) exp {-iojzjc) is a slowly varying function of z on the scale 
of the wavelength A, and for the sake of simplicity will also be named "the 

field". In Eq. (17) of reference [ITfl an expression for E{z,r,iXi) generated by 
an electron moving along a generic trajectory r(z) was found. Working out 
that equation under the resonance approximation for the case of a planar 

undulator where r(z) = K/{yk^) cos (fcjyz)^^ + I + zrf, Cx being the unit vector 
in the x-direction, yields the horizontally polarized field 



{-e)Ka)Ajj 
2c2y 



L„,/2 

r dz' 

J z-z' 



exp 



-L„,/2 



c + 



2c 



z' + 



2c{z - z') 



(22) 



Here we defined the detuning parameter C = oj/{2ylc) - k-a, = {Aa)/ixir)kju, 
where co = co,. + Aco. Thus, C specifies "how much" co differs from the 
fundamental resonance frequency a)r = 2y^c/Aa,. The subscript "C" in 
Ec indicates that this expression is valid for arbitrary detuning param- 
eter. Moreover, Az„ = l/k^u = KoK^tl), Xj„ being the undulator period; 

= y/{l + K^/2); K = {Ay^eH^a) I {2nmeC^) is the undulator parameter, 
being the electron mass and being the maximum of the magnetic field 
produced by the undulator on the z axis. Finally, Lu, is the undulator length 
andAjj = /o[KV(4 + 2^2)] -/^[KV (4 + 2K2)]j^ indicating the Bessel function 
of the first kind of order n. It should be stressed that Eq. (23) was derived un- 
der the resonance approximation meaning that the large parameter N^, » 1 
was exploited, together with conditions Aojjojr <?c 1 andC-i-(i'r^/(2cz^) <?c k-,„, 
meaning that we are looking at frequencies near the fundamental and an- 
gles within the main lobe of the directivity diagram. Moreover, the reader 
should keep in mind that no focusing elements are accounted for in the 
undulator. This fact is intrinsically related to the choice of r(z) done above. 



Further algebraic manipulations (see Appendix B of |130) show that Eq. (22) 
can be rewritten as: 
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L„/2 

~ {-e)Ka)Ajj f dz' 
£c = 



2c^y J z-z' 

-L„,/2 



X exp 



2cz 



2c(z - z') 



z z 



(23) 



In this paper we will make a consistent use of dimensional analysis, which 
allows one to classify the grouping of dimensional variables in a way that 
is most suitable for subsequent study. Normalized units will be defined as 



E = -- — ^—E, 



KcoeA 



JJ 



n = n 



IojL,, 



C = Lu,C = 2nN, 



0) - COr 



OJr 



r = r 



l = L 



0) 



L c ' 



CO 



L c ' 



1-171 



(24) 



Moreover, for any distance z, we introduce 6 = f/z. The algorithm for calcu- 
lating the cross-spectral density will be formulated in terms of dimensionless 
fields. Therefore we re-write Eq. (23) as 



Ec = exp 
where 



z 
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z, C, 



0---T7 
z ' 



(25) 



1/2 

f dz' 

Wc{z,C,a)= I T^Tjexp 



Cz' + 



zz a"^ 
2(z - z') 



(26) 



-1/2 



A physical picture of the evolution of the field along the z direction was 
given in reference [17], where Fourier optics ideas were used to develop a 
formalism ideally suited for the analysis of any SR problem. In that reference. 
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the use of Fourier optics led to establish basic foundations for the treatment 
of SR fields, and in particular of undulator radiation, in terms of laser 
beam optics. Radiation from an ultra-relativistic electron can be interpreted 
as radiation from a virtual source, which produces a laser-like beam. In 
principle, such virtual source can be positioned everywhere down the beam, 
but there is a particular position where it is similar, in many aspects, to the 
waist of a laser beam. In the case of an undulator this location is the center of 
the insertion device. A virtual source located at that position ("the" virtual 
source) exhibits a plane wavefront. Therefore, it is completely specified by 
a real-valued amplitude distribution of the field (see Eq. (34) of Jl^). This 
amplitude can be derived from the far zone field distribution. Free-space 
propagation from the virtual source through the near zone and up to the 
far-zone, can be performed with the help of the Fresnel formula: 

E{z,r) = TTl^—TT I df' E{Zs,f')exp 
2n{z - Zs) J 



r' 



2{z - Zs) 



(27) 



where the integral is performed over the transverse plane, and Zg is the 
virtual source position down the beamline. 

These considerations were applied in [17] to the case of undulator radiation 
under the applicability region of the resonant approximation. With reference 
to Fig. 6, we let z = be the center of the undulator. Thus, the position of the 
virtual source is fixed in the center of the undulator too, Zg = 0. For simplic- 
ity, the resonance condition with the fundamental harmonic was assumed 
satisfied, i.e. a; = cu,-. For this case, an analytical description of undulator 
radiation was provided. The horizontally-polarized field produced by a sin- 

gle electron with offset I and deflection r) in the far-zone (i.e. at z » 1) can 
be represented by the scalar quantit}© 



I z, //, /, 1 = I exp 



./z02 ^ ^ p 

i\ 6-1 + — - 

\ 2 2z2 



(28) 



^ Note that for a particle moving on axis, at / = and ?] = the quadratic phase 
exp[/z0^/2] in Eq. (28) is indicative of a spherical wavefront in paraxial approxima- 

tion on the observation plane. When / is different from zero, the laser-like beam is 

shifted, and this justifies the present of extra-factors including Z. When the particle 
also has a deflection fj, the laser-like beam is tilted, but the wavefront remains 
spherical. Since the observation plane remains orthogonal to the z axis, the phase 
factor before Wf does not include t] and thus it does not depend on the combination 

S - II z - fi- 
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Fig. 6. Illustration of the undulator geometry and of the observation plane, 
where 



^f{a) = sine 



T 



(29) 



subscript / indicating the "far-zone". The field distribution of the virtual 
source positioned at z = 0, corresponding to the waist of our laser-like beam 
was found to be: 



where 



-in exp 



f-l 



(30) 



%(a) = ^[7i-2Si(a2)] . (31) 

where Si(z) = dtsm{t)/t indicates the sin integral function and subscript 
"0" is indicative of the source position. Plots of and ^re given in 
Fig. 7 and Fig. 8. It should be noted here that the independent variable in 
both plots is the dummy variable a. The characteristic transverse range of 
the field in the far zone is in units of the radiation diffraction angle yfKjL^, 
while the characteristic transverse range of the field at the source is in units 
of the radiation diffraction size yfKL^. 

Finally, with the help of the Fresnel propagation formula Eq. (27), we found 
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Fig. 7. Universal function ^/(a), used to calculate the far-zone radiation field of a 
single electron at the fundamental harmonic at perfect resonance. 
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Fig. 8. Universal function ^o(^^)/ used to calculate the radiation field of a single 
electron at the fundamental harmonic at perfect resonance on the virtual source 
plane. 
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the following expression for the field distribution at any distance z > 1/2 
from the virtual source: 



exp 



z 
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z, 



e---r, 



(32) 



where we defined 



W{z, a) = exp 



.za"^ 




' ■'^2 2 " 

iz^a^ 




' ■'^2 2 " 

iz^a^ 






-Ei 


1 




2z- 1 


2z + l 





(33) 



Xoo 
^ tif exp(-f ) / 1 indicates the exponential integral function. Eq. 

(32) is a particular case of Eq. (25) at perfect resonance. Note that free space 

basically acts as a spatial Fourier transformation. This means that the field 

in the far zone is, aside for a phase factor, the Fourier transform of the field 

at any position z down the beamline. It is also, aside for a phase factor, the 

spatial Fourier transform of the virtual source: 



r "* 1 

^ I iz r z|0p ^ I ^ ^\ 

£o|0,r/,/,f 1 = -— I d0exp — —z Ef\z,f\,l,e\exY> 



if • 6 



(34) 



It follows that 



00 

^ dd ^/(0) exp [zf ■d\=ln J dd djo (\f\d) ^/(0) = 271^^0 (l^l) • (35) 



We conclude verifying that Eq. (32) is in agreement with Eq. (28) and Eq. (30) 
z » 1 and for z = 0, respectively. Consider first z - 0. For positive numbers 

2 



r-l-ZT] 



> 0, we have 



- i[n + 2Si(a^)] = Ei(-za^) - Ei(fa^) , 
and Eq. (32) yields back Eq. (30). 



(36) 



Consider now the limit for z » 1. For positive numbers 
one has 



d-fj 



> 



exp 



zz^^ 



2a2 



Ei 



' iez 

.2z-l. 
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Ei 



,2z + l 



isinc(f), (37) 
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as it can be directly seen comparing Eq. (32) with Eq. (25)), where integration 
in Eq. (26) is performed directly at C = and in the limit for z » 1. Thus, 
Eq. (32) yields back Eq. (28), as it must be. 

Finally, it should be noted that expressions in the present Section 2.4 have 
been derived for co > 0. Expressions for the field at negative values of co 

can be obtained based on the property E{-aji) = E {co) starting from explicit 
expressions at a; > 0. 

2.5 Cross-spectral density of an undulator source and its free-space propagation. 



From now on we consider a normalized expression of the cross-spectral 
density, G, that is linked with G in Eq. (12) by a proportionality factor 



Moreover, we introduce variables 

-> -> 

. ^ ^ ^ ,r.c,\ 

Ar - r-i-r2 , r - — - — . (39) 



and 



Ad = 61-62, 6= \ \ (40) 



where, as before, 6 = f/z. Thus, Eq. (12) can now be written as 

G(z,f,Ar:C) = (£ /],/,z,f + y,C F /),/,z,f - y,cU . (41) 



On the one hand, the cross-spectral density as is defined in Eq. (41) includes 
the product of fields which obey the free space propagation relation Eq. (27). 
On the other hand, the averaging over random variables commutes with 
all operations involved in the calculation of the field propagation. More 
explicitly, introducing the notation £(z) = (9[£(Zs)] as a shortcut for Eq. (27) 
one can write 



G(z) = {O [£(z,)] O* [e\Zs)\i =0 0* [(e(z,)F(z,)}] =0a [G(z,)] , (42) 
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where O may also represent, more in general, any linear operator. Once the 
cross-spectral density at the source is known, Eq. (42) provides an algorithm 
to calculate the cross-spectral density at any position z down the beamline 
(in the free-space case). Similarly, propagation through a complex optical 
system can be performed starting from the knowledge of G(Zs). As a result, 
the main problem to solve in order to characterize the cross-spectral density 
at the specimen position is to calculate the cross-spectral density at the 
virtual source. For the undulator case, we fix the position of the source in 
the center of the undulator Zg = 0. This is the main issue this paper is devoted 
to. However, free-space propagation is also treated, and may be considered 
an illustration of how our main result can be used in a specific case. 

Based on Eq. (41) and on results in Section 2.4 we are now ready to present an 
expression for the cross-spectral density at any position down the beamline, 
always keeping in mind that the main result we are looking for is the cross- 
spectral density at the virtual source position. 

We begin giving a closed expression for C valid at any value of the detuning 
parameter C by substituting Eq. (25) in Eq. (41), and replacing the ensemble 
average with integration over the transverse beam distribution function. We 
thus obtain 



X 



|z,C,0,A0j = 



exp 



izQ ■ A0 



d\ exp 



-il ■ Ad 



Z,Cr 



2 z ' 









z, C, 







^ Ad I ^ 
9 - -r- - - - ri 
2 z ' 



(43) 



It is often useful to substitute the integration variables / with (p = -O + l/z + i] 
In fact, in this way, Eq. (43) becomes 



'^2 

z exp 



;(z,C,0,A0j = 



izd • Ad 



J' dcp J' df] 



exp 



z, C, 



Ad 



iz^^ + d-^yAd 
\ 

(44) 



Ad 



For choices of / of particular interest (e.g. product of Gaussian functions for 
both transverse and angle distributions), integrals in dff can be performed 
analytically, leaving an expression involving two integrations only, and still 
quite generic. 



Eq. (44) is as far as we can get with this level of generality, and can be 
exploited with the help of numerical integration techniques. However, it 
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still depends on six parameters at least: four parameter£3 are needed to 
specify fj_, plus the detuning parameter C and the distance z. 

In the following we will assume C <^ 1, that allows us to take advantage 
of analytical presentations for the single-particle field obtained in [flT}] and 
reported before. This means that monochromatization is good enough to 
neglect finite bandwidth of the radiation around the fundamental frequency 
By this, we automatically assume that monochromatization is performed 
around the fundamental frequency. It should be noted, however, that our 
theory can be applied to the case monochromatization is performed at other 
frequencies too. Analytical presentation of the single-particle field cannot be 
used in full generality, but for any fixed value C of interest one may tabulate 
the special function once and for all, and use it in place of W throughout 
the papeiP^ From this viewpoint, although the case of prefect resonance 
studied here is of practical importance in many situations, it should be 
considered as a particular illustration of our theory only. 

Also note that in Eq. (43) the electron beam energy spread is assumed to be 
negligible. Contrarily to the monochromator bandwidth, the energy spread 
is fixed for a given facility: its presence constitutes a fundamental effect. 
In order to quantitatively account for it, one should sum the dimensionless 
energy-spread parameter Ag = AnN-,v6y/y to C in Eq. (43) and, subsequently, 
integration should be extended over the energy-spread distribution. Typical 
energy spread 5y/y for third generation light sources is of order 0.1%. For 
ERL sources this figure is about an order of magnitude smaller, 6y/y ~ 
0.01%. 

In order to study the impact of a finite energy spread parameter Ag, of 
a finite radiation bandwidth and of a relatively small detuning from the 
fundamental we consider an expression for the intensity of a diffraction- 
limited beam z^I = z^{|Ep) including both Ag and C: 



z2j(0, CAe) = 



V27iAi 



oo 

/ 



d^E exp 



2A^ 



(45) 



We plotted z^J for different values of C and Ae in Fig. 9. First, in Fig. 9 
(a), we compared, at C = 0, the case for negligible energy spread with 



^ At least. This depends on the number of parameters needed to specify f±. For a 
Gaussian distribution in phase space, four parameters are needed, specifying rms 
transverse size and angular divergence in the horizontal and vertical direction. 

Of course, selection of a particular value C still implies a narrow monochromator 
bandwidth around that value. 
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Fig. 9. Study of I given in Eq. (45) for different valued of C and Af. Here N-a, - 70. 
Plot (a) : comparison, at C = 0, of the case for negligible energy spread with the 
case A£ = 0.88. Plot (b) : comparison, at negligible energy spread, of the case C = 
with the case C = ±0.5. Plot (c) : comparison of the case C = at negligible energy 
spread with the case C = +0.5 at Aco/co ^ 0.1%. 

the case = 0.88, corresponding to 6y/y ^ 0.1% at Nj„ = 70, typical 
of third generation sources. As one can see, maximal intensity differences 
are within 10%. Second, in Fig. 9 (b) we compared, at negligible energy 
spread, the case C = with C = ±0.5 at N^, = 70, corresponding to a shift 
Aco/co ^ 0.1%. One can see that also in this case maximal intensities differ of 
about 10%. Analysis of Fig. 9 (c), where the case C = at negligible energy 
spread is compared with cases C = ±0.5 at Aco/co ^ 0.1% and Ni„ = 70 
leads to a similar result. This reasoning allows to conclude that our simplest 
analytical illustrations can be applied to practical cases of interest involving 
third generation sources and undulators with up to 70 periods with good 
accuracy. It should be remarked that such illustration holds for the first 
harmonic only. In fact, while the shape of z^I is still given by Eq. (45) in 
the case of odd harmonic of order h, parameters C and Ae are modified 
according to Ch - hC and Aeh = hA^, decreasing the applicability of our 
analytical results. 

With this in mind, we can present an expression for the cross-spectral density 
at C <^c 1 based on Eq. (41) and Eq. (32). Substituting the latter in the former 

we obtain an equation for G(z, 6, Ad) that can be formally derived from 
Eq. (43) by substitution of Wc with W. Similarly as before, one may give 

alternative presentation of G(z, 6, Ad) replacing the integration variables I 

^ -t -> ^ ^ ^ 

with cp = -6 + 1/2 + 1]. This results in another expression for G(z, 6, Ad) that 

can be formally derived from Eq. (44) by substituting Wc with W. This last 
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expression presents the cross-spectral density in terms of a convolution of 
the transverse electron beam phase space distribution with an analytical 
function, followed by Fourier transformatiorl"!. 

One may obtain an expression for G at z = as a limiting case of Eq. (43) or 
Eq. (44) at C = 0. It is however simpler to do so by substituting Eq. (30) in 
Eq. (41) that gives 



G ^0, r, Ar^ = J' dif exp [Of ■ Ar\ 



2 



r + I 

2 



(46) 



where the function has already been defined in Eq. (31). The product 
^0 (|^~ Ar/2|^ ^0 (|^+ Ar/2|^ is a four-dimensional , analytical function in 
f and Ar. Eq. (46) tells that the cross-spectral density at the virtual source 
position can be obtained convolving (jf- Ar/2|^ (|^+ Ar/2|^ with the 

transverse beam distribution at z = 0, /j_ Z, r) , considered as a function of I, 



and taking Fourier transform with respect to f]. When the betatron functions 
have minima in the center of the undulator we have 



Then, Eq. (46) becomes 



(47) 



G(0,r,Ar^ = ^ d^/,, (^) exp [z^ • Ar] 



X 



Af_. 
2 



1% 



^ Ar i 

r+ / 

2 



(48) 



that will be useful later on. In this case, the cross-spectral density is the 
product of two separate factors. First, the Fourier transform of the distri- 
bution of angular divergence of electrons. Second, the convolution of the 
transverse electron beam distribution with the four-dimensional function 
vFo(|r-Ar72|)vi/o(|r + Ar72|). 

Eq. (43) (or Eq. (44)) constitutes the most general result in in the calculation 
of the cross-spectral density for undulator sources. Its applicability is not 



Aside for an inessential multiplicative constant. This remark also applies in what 
follows. 
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restricted to third generation light sources. In particular, it can be used for 
arbitrary undulator sources like ERLs [1] or XFEL spontaneous undulators 
ifisll . Eq. (44) has been derived, in fact, under the only constraints » 1, 
» 1 and ajCOr » Njy. Note that ctj ~ 30 ps for a typical SR source, 
whereas Oj ~ 100 fs for an XFEL spontaneous undulator source or an ERL. 
Yet, for all practical cases of interest, OjOJ » N^. As we have seen before, 
Eq. (44) further simplifies in the particular but practical case of perfect 
resonance, i.e. in the limit for C «; 1. A particularly important asymptote of 
Eq. (43) at perfect resonance is at the virtual source position, described by 
Eq. (46), which express the cross-spectral density in the undulator center. 
While Eq. (43) (or Eq. (44)) solves all problems concerning characterization of 
transverse coherence properties of light in free-space, the knowledge of Eq. 
(46) constitutes, in the presence of optical elements, the first (and main) step 
towards the characterization of SR light properties at the specimen position. 
In fact, the tracking of the cross-spectral density can be performed with the 
help of standard statistical optics formalism developed for the solution of 
problems dealing with partially coherent sources. Finally, it should be noted 
that in the case of XFELs and ERLs, there is no further simplification that 
we may apply to previously found equations. In particular, the transverse 
electron beam phase space should be considered as the result of a start- 
to-end simulation or, better, of experimental diagnostics measurements in 
a operating machine. On the contrary, as we will see, extra-simplifications 
can be exploited in the case of third-generation light sources, allowing for 
the development of a more specialized theory. 

Inspection of Eq. (43) or Eq. (46) results in the conclusion that a Gaussian- 
Schell model cannot be applied to describe partially coherent SR light. In 
fact, functions Wq, M^o and are of non-Gaussian nature, as the laser-like 
beam they can be ascribed to is non-Gaussian. This explains our words in the 
Introduction, where we stated that are of general theoretical interest, 

but they do not provide a satisfactory approximation to third generation SR 
sources. 

As a final remark to this Section, we should discuss the relation of our 
approach with that given, in terms of Wigner distribution, in [[ill, 12]. As 



said in Section 1, treatment based on Wigner distribution is equivalent to 
treatment based on cross-spectral density. We chose to use cross-spectral 
density because such quantity is straightforwardly physically measurable, 
being related to the outcome of a Young's experiment. Essentially, one can 
obtain a Wigner distribution W from G by means of an inverse Fourier 
transformation: 

W (z, ^'^) = ^ J d{Ar) exp [-iAr- u] G (z, r, Ar) . (49) 
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Thus, Eq. (48) gives 



w(o,f,fi) = jdtj 4/(f,^)w„(o, 



r -l,u-r\ 



(50) 



where 



Wo (O, a, 6) = i rf(Ar ) exp [-i^r■ b\ 



a - 



« H 

2 



, (51) 



The Wigner distribution at z = is presented as a convolution product 
between the electron phase-space and a universal function Wp. This result 



may be directly compared (aside for different notation) with [111, 11211 , where 



the Wigner distribution is presented as a convolution between the electrons 
phase space and a universal function as well. The study in |lll, Il2|1 ends at 
this point, presenting expressions for arbitrary detuning parameter. On the 
contrary, in the following Sections we will take advantage of expressions at 
perfect resonance, of small and large parameters related to third generation 
light sources and of specific characteristics of the electron beam distribution. 
This will allow us to develop a comprehensive theory of third generation 
SR sources. 



3 Theory of transverse coherence for third-generation light sources 



3.1 Cross-Spectral Density 



We now specialize our discussion to third-generation light sources. We as- 
sume that the motion of particles in the horizontal and vertical directions 
are completely uncoupled. Additionally, we assume a Gaussian distribu- 
tion of the electron beam in the phase space. These two assumptions are 
practically realized, with good accuracy, in storage rings. For simplicity, we 
also assume that the minimal values of the beta-functions in horizontal and 
vertical directions are located at the virtual source position z = 0, that is 

often (but not alwayQ the case in practice. /j_ = f,^SVx)UiVy)hi^x)Ai^y) 
with 



exp 



'2D, 



^2^ 



: exp 



'2D. 



Generalization to the case when this assumption fails is straightforward. 
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Here 



n = z_ 



(53) 



cr;c,y and Ox',y' being rms transverse bunch dimensions and angular spreads. 
Parameters Nx,y will be indicated as the beam diffraction parameters, are 
analogous to Fresnel numbers and correspond to the normalized square of 
the electron beam sizes, whereas D^^y represent the normalized square of the 
electron beam divergences. Consider the reduced emittances ex,y - ex,y/A, 
where ex,y indicate the geometrical emittance of the electron beam in the 
horizontal and vertical directions. Since we restricted our model to third 
generation light sources, we can consider » 1. Moreover, since betatron 
functions are of order of the undulator length, we can also separately accept 



» 1 , Dx » 1, (54) 

still retaining full generality concerning values of Ny and Dy, due to the 
small coupling coefficient between horizontal and vertical emittance. 

Exploitation of the extra-parameter e^; » 1 (or equivalently » 1 and 
Dx » 1) specializes our theory to the case of third-generation sources. 

With this in mind we start to specialize our theory beginning with the 
expression for the cross-spectral density at the virtual source, i.e. Eq. (48). 

After the change of variables (p — > -r + 1, and making use of Eq. (52), Eq. 
(48) becomes 
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where Ax and Ay indicate components of Ar. In Eq. (55) the range of variable 
(px is effectively limited up to values \(px\ ~ 1. In fact, (px enters the expression 
for It follows that at values larger than unity the integrand in Eq. (55) 
is suppressed. Then, since Nx » 1, we can neglect (px in the exponential 
function. Moreover Dx ^ 1 and from the exponential function in follows 
that Ax « 1 can be neglected in M^o- As a result, Eq. (55) is factorized in 
the product of a horizontal cross-spectral density Gx{0, x. Ax) and a vertical 
cross-spectral density Gy(0, x. Ax): 



G (O, r, Arj = Gx (0, x. Ax) Gy (0, y. Ay) , 



(56) 



where 
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(57) 
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(58) 



Note that, in virtue of Eq. (27), factorization holds in g;eneral, at any posi- 
tion z. This allows us to separately study Gx and Gy. Gx describes a quasi- 
homogeneous Gaussian source, which will be treated in Section 4.1. Here 
we will focus our attention on Gy only. It should be remarked that the quasi- 
homogenous Gaussian source asymptote is obtained from Eq. (58) in the 
limit » 1 and Dy » 1. In other words, normalization constants in Eq. 
(57) and Eq. (58) are chosen in such a way that Eq. (58) reduces to Eq. (57) 
in the limit Ny » 1 and Dy » 1 (with the obvious substitution x — > y). It 
should be clear that this normalization is most natural, but not unique. The 
only physical constraint that normalization of Eq. (57) and Eq. (58) should 
obey is that the product GxGy should not change. 



Let us define the two-dimensional universal function S{a, 6) aa 



S stands for "Source'' 
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Fig. 10. Plot of the universal function S, used to calculate the cross-spectral density 
at the virtual source when Nx ^ I, l,Ny and Dy are arbitrary. 



(59) 



The normalization constant TCg ^ 0.714 is chosen in such a way that 3(0, 0) = 
1. We can present Cy at the virtual source with the help of S as 







oo 

J d(py exp 

— OO 


((pv + y) 
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2Ny 



x^^'S{cPy,Ay) 



(60) 



Therefore, Gy at the virtual source is found by convolving an universal 
function, S, with a Gaussian function and multiplying the result by another 
Gaussian function. 



Note that in the limit N;^ » 1 and D;^ » 1 there is no influence of the 
electron beam distribution along the vertical direction on Gx. In spite of 
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Fig. 11. Plot of the universal function , used to calculate the cross-spectral density 
in the far zone when tiy^ » 1, D;^ » 1, Ny and Dy are arbitrary. 



this, in the same limit, Eq. (60) shows that there is an influence of the 

horizontal electron beam distribution on Gy, due to the non-separability of 
the function M^o in S. In fact, contrarily to the case of a Gaussian laser beam, 
%[(/);,^ -f (a - 5/2)2] ^ m^\(^^Y^^\{a - 5/2)^], and the integral in d(^y„ that is 
a remainder of the integration along the x-direction, is still present in the 
definition of S. However, such influence of the horizontal electron beam 
distribution is independent of N;c and D^. As a consequence, tS is a universal 
function. 

A plot of S is presented in Fig. 10. <S is a real function. Then, G at the virtual 
plane is also real. Moreover, S is invariant for exchange of a with 6/2. 

Let us now deal with the evolution of the cross-spectral density Gy in free- 
space. In principle, one may use Eq. (60) and apply Eq. (42), remembering 
Eq. (27). It is however straightforward to use directly Eq. (44) at C «: 1. 
Under the assumption N;^ ^ 1 arid D;^ » 1, as has been already remarked, 
factorization of (5 as a product of 0,^ and (5y holds for any value of z. Isolating 
these factors in Eq. (44) at C «; 1 and using Eq. (52) one obtains 
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Gy (Z, By, Ady) = CXp [zZSyAej ^ d(/)y ["ZZ + 0^ AO y] 



X 
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J" ir)y exp |^zz/)yA0yj exp 




2Ny/z2 



(Pl + \(py + 



l/2^ 



• (61) 



The integral in dfjy (second line) can be calculated analytically yielding 




2nypR^Ny/z^ + Dy 

i{c(>y + dy)AeyNy/Z 



Ny/Z^+Dy 



exp [zf0yA0yj exp 



2{Ny/Z^+Dy) 



exp 



2(Ny/z2+Dy) 



(pi + \(py- 



A0. 



(t>l + \(py + 



A0„ 



l/2^ 



(62) 



We are now interested in discussing the far-zone limit of Eq. (62). Up to 
now we dealt with the far-zone region of the field from a single particle, Eq. 
(28). In this case, the field exhibits a spherical wavefront. Such wavefront 
corresponds to the quadratic phase factor in Eq. (28). Note that when the 
electron is moving on-axis, Eq. (28) consists of the product of z"^ exp[fz0^/2] 
by a real function independent of z. Such field structure can be taken as 
a definition of far-zone. A similar definition can be used for the far-zone 
pertaining the cross-spectral density. We regard the quadratic phase factor 
exp[zz0yA0y] in Eq. (62) as the equivalent, in terms of cross-spectral density, 
of the quadratic phase factor for the single-particle field. We therefore take 
as definition of far-zone the region of parameters where Gy(z, By, AOy) = 
z~^h{By, Ady) exp[zz0yA0y], k being a real function, and remains like that for 
larger values of z. 

Let us discuss the definition of the far-zone region in terms of problem pa- 
rameters. In the single-particle situation, the only parameter of the problem 
was z and, as is intuitively sound, the far-zone region was shown to coincide 
with the limit z » 1. In the case of Eq. (62), we deal with three parameters 
z, Ny and Dy. Therefore we should find that the far-zone is defined in terms 
of conditions involving all three parameters. 

When Ny < 1 and Dy < 1, analysis of Eq. (62) shows that the far-zone region 
is for z » 1. In this case the definition of far-zone for Gy coincides with that 
of far-zone for the field of a single particle. 
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However, when either or both Ny » 1 or Dj^ » 1 the situation is different, 
and one finds that the far-zone condition is a combination of z, Ny, and 
Dy. In all these cases, detailed mathematical analysis of Eq. (62) shows that 
far-zone is reached when 

max[Ny, 1] « max[Dy, 1] . (63) 

As it will be clearer after reading Section 4.4, but can also be seen considering 
the definition of our dimensionless units, the physical meaning of compar- 
isons of Ny and Dy with unity in condition (63) is that of a comparison be- 
tween diffraction-related parameters (diffraction angle and diffraction size) 
and beam-related parameters (divergence and size of the electron beam). 

When Ny » 1, but Dy <1 condition (63) reads » » 1. This result is in 
agreement with intuition. The far-zone condition for Gy does not coincide 
with that for the field of a single electron, but it is anyway reached far away 
from the source, at z » 1. As we will see in Section 4.2, the case Ny » 1 
with Dy <1 corresponds to a quasi-homogeneous non-Gaussian source. In 
Section 4.2 we will see that in this case the VCZ theorem is applicable, and 
its region of applicability is in agreement with our definition of far-zone 

Z2 » Ny » 1. 

When Dy » 1 (with arbitrary Ny) and condition (63) holds, analysis shows 
that the phase factor under the integral sign in Eq. (62) can be neglected. 
Moreover, in this case, the Gaussian function in (py + dy has a width in (py 
much larger than unity, while the integral in dcp^ in second line of Eq. (62) 
has a width in (py of order unity, because does not depend on parameters. 
Therefore, the dependence in c^y in the Gaussian function can always be 
omitted, and the Gaussian function factors out of the integral sign in dcpy. 
One is left with the product of exponential functions and the double integral 



/(A0y) ^ 




which has a very peculiar property. In fact, it does not depend on z . The proof 
is based on the autocorrelation theorem in two dimensions, and is given in 
detail in Appendix C of reference [13]. This quite remarkable property of / 
carries the consequence that substitution of with can be performed in 
Eq. (64) without altering the final result. 



36 



When both Dy » 1 and Ny » 1, one may neglect the dependence in Ady 
in functions W and in Eq. (64), because the exponential function in A0y 
before the integral sign limits the range of Ady to values of order 1 / ^jNy ^ 1 . 
As a result, the double integration in dcpx and dcpy yields a constant, and the 
description of the photon beam is independent of W, i.e. does not include 
diffraction effects. This result is intuitive: when the electron beam size and 
divergence is large compared to the diffraction size and divergence, the 
photon beam can be described in terms of the phase-space distribution of 
the electron beam. This approach will be treated in more detail in Section 
4.4. 

Finally, when Dy » 1 and Ny < 1, diffraction effects cannot be neglected, 
nor can be the dependence in A0y in W and M^*. In this case, from condition 
(63) we obtain that the far-zone coincides with DyZ^ » 1. This result is 
completely counterintuitive. In fact, since Dy » 1, the far-zone is reached 
for values z ~ 1, i.e. at the very end of the undulator. Yet, diffraction effects 
cannot be neglected, and the field from a single electron is far from exhibiting 
a spherical wavef ront at z ~ 1 . This paradox is solved by the special property 
of the double integral in Eq. (64), that allows one to substitute with Wf 
independently of the value of z. 

As it will be discussed in Section 4.1 and Section 4.2, the case Dy » 1 
corresponds to a quasi-homogeneous Gaussian source when Ny » 1 and to 
a quasi-homogeneous non-Gaussian source when < 1. It will be shown 
that the VCZ theorem is applicable to these situations. In particular, the 
applicability region of the VCZ theorem will be seen to be in agreement 
with our definition of far-zone. 

Our discussion can be summarized in a single statement. The far zone is 
defined, in terms of problem parameters, by condition (63). Remembering 
this condition one can derive the far-zone expression for Gy simplifying Eq. 
(62): 



1 1 

Gy(z, dy, Ady) = =exp [K0yA0y] exp 

z 2n ^JnDy 



NyAdl 



X J #yexp -^^^^^ ^/[^y^^y) ' (65) 

— CO 

where the universal two-dimensional functiorQ3 'J^{oc, 5) is normalized in 
such a way that fiO, 0) = 1 and reads : 



14 r stands for "Far". 
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(66) 



withTCf = 3/(8 Vti). 



Note that, similarly to the source case, the right hand side of Eq. (65) is 
found by convolving an universal function, 9^, with a Gaussian function 
and multiplying the result by another Gaussian function. 

A plot of ;F function defined is given in Fig. 11. !F is a real function. Thus, 
only the geometrical phase factor in Eq. (65) prevents Gy in the far-zone from 
being real. Another remarkable property of is its invariance for exchange 
of a with 6/2. Also, is invariant for exchange of a with -a (or 6 with -6). 

Finally, from Eq. (35) and Eq. (59) we have 



S{x, y) 



oo oo 

J' da J' d6 exp 
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(67) 



3.2 Intensity distribution 



With expressions for the cross-spectral density at hand, it is now possible 
to investigate the intensity distributiorQll along the beamline, letting Ar = 
in the expression for G. Since Nx » 1 and » 1, factorization of the 
cross-spectral density still holds. Therefore we can investigate the intensity 
profile along the vertical direction without loss of generality. 

Posing Ay = in Eq. (60) we obtain the intensity profile at the virtual source, 
ly (0, y), as a function of y: 



4(0'y) = 




yexp 



(^y + yf 



2N. 



(68) 



15 Words "intensity distribution" include some abuse of language here and in the 
following. What we really calculate is the ensemble average of the square modulus 
of the normalized field, Ixly = Conversion to dimensional units, followed 

by multiplication by c/(47z^) yields the spectral density normalized to the electron 
number Ng. 
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Fig. 12. The universal function Ig, used to calculate intensity at the virtual-source 
position. 

where we introduced the universal function 



00 

Isia) = S{a, 0) = %s j #x^o ( + • (69) 

—oo 

A change of the integration variable: (px — > x = {(pi + a^)^^^ allows the 
alternative representation: 

oo 

Isia) = TKs r dx'-^^^M^^l (x) , (70) 



where the function rect(^) is defined following |[^: it is equal to unity for 
1^1 ^1/2 and zero otherwise. 

The intensity at the virtual-source position is given in terms of a convolution 
of a Gaussian function with the universal function Is- A plot of Is is given 
in Fig. 12 . 
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Fig. 13. The universal function Jf , used to calculate the intensity in the far zone. 

Similar derivations can be performed in the far zone. Posing Kdy = in Eq. 
(65) we obtain the directivity diagram of the radiation as a function of By. 



oo 



2D, 



2-| 



(71) 



where we defined 



1/2 



If {a) = r {a, 0) = TCf J {(^l + a') 

— oo 
oo 



(72) 



The intensity in the far-zone is given in terms of a convolution of a Gaussian 
function with the universal function Jf. A plot of If is given in Fig. 13. 
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3.3 Spectral degree of coherence 



We can now present expressions for the spectral degree of coherence at the 
virtual source (that is a real quantity) and in the far-zone (that is not real). 
Eq. (14) can be written for the vertical direction and in normalized units as: 



Gy{z,y,Ky) 



[ly{z, y + Ay/2)f' [lyiz, y - Ay/2)] 



1/2 • 



(73) 



Substitution of Eq. (60) and Eq. (68) in Eq. (73) gives the spectral degree of 
coherence at the virtual source: 



gy(0,y,Ay) = exp\- 
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(74) 



Similarly, substitution of Eq. (65) and Eq. (71) in Eq. (73) gives the spectral 
degree of coherence in the far zone: 



gy{Z, dy, Ady) = eXp [zZ0yA0y] eXp 
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(75) 
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3.4 Influence of horizontal emittance on vertical coherence for Dy <^ 1 and 



Ny«l 



The theory developed up to now is valid for arbitrary values of Ny and 
Dy. In the present Section 3.4 we discuss an application in the limiting case 
for Dy <^ 1 and Ny <^ 1 corresponding to third generation light sources 
operating in the soft X-ray range. 

At the virtual source position, the following simplified expression for the 
spectral degree of coherence in the vertical direction is derived from Eq. 
(74): 

gy(0,y,Ay)^Xs{y,Ay), (76) 
where 

S{a, 6) 

[Is{a - 6/2)f'^[Is{a + 6/2)]' 



Thus, gy is given by the universal function Xs- A plot of Xs is given in Fig. 
14. 

Similarly, in the far zone, one obtains from Eq. (75): 



gy{Z, dy, Ady) = eXp [izd y Ad y\X ^{6 y , Ad y) . (78) 

Here the universal function Xf is given by 

^Fia, 6) = ^ . (79) 

Thus, \gy\ is given by the universal function Xf- A plot of Xf is presented in 
Fig. 15 as a function of dummy variable a and 6. It is straightforward to see 
that XF{dy, Ady) is symmetric with respect to Ady and with respect to the 
exchange ofA0y/2 with 0y. When 0y = 0,i.e. dyi - -0y2/ we obviously obtain 
/Yf(0, A0y) = 1 that corresponds to complete coherence at this particular 
value of By. However, since = Xf oscillates from positive to negative 
values, in general one never has full coherence in the vertical direction, even 
in the case of zero vertical emittance. Note that this effect does not depend, 
in the limit for Nx^l and Dx » 1, on the actual values of Nx and Dx. 

In other words, as it is evident from Fig. 15, XFiOy,Ady) exhibits many 
different zeros in 6y for any fixed value of Ady. In Fig. 16 some of these 
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Fig. 14. Plot of the universal function Xs, used to calculate the modulus of the 
spectral degree of coherence of the source when Nj » 1, Dj » 1, Ny <sc 1 and 

Dy « 1. 

zeros are illustrated with black circles on the plane {By, Ady). Consider a 
two-pinhole experiment as in Fig. 4. Once a certain distance zAdy between 
the two pinholes is fixed. Fig. 16 illustrates at what position of the pinhole 
system. By, the spectral degree of coherence in the vertical direction drops 
from unity to zero for the first time. 

To estimate the importance of this effect, it is crucial to consider the position 
of By in the directivity diagram of the radiant intensity, that coincides in 
this cas with IpiBy) (solid line in Fig. 16). From Fig. 16 one can see that 
XpiBy, ABy) drops to zero for the first time at ABy ~ 2 0j/ ~ 2, where the 
X-ray flux is still intense. This behavior of the degree of coherence may 
influence particular kind of experiments. To give an example we go back 
to the two-pinhole setup in Fig. 4. After spatial filtering in the horizontal 
direction, one will find that for some vertical position By of the pinholes 
(at fixed ABy), well within the radiation pattern diagram, there will be no 
fringes, while for some other vertical position there will be perfect visibility. 



° In other words, it can be shown that Jf is the directivity diagram corresponding 
to the case » 1, D;^ » 1, Ny <sc 1, Dy «c 1. 
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Fig. 15. Plot of the universal function X^, used to calculate the modulus of the 
spectral degree of coherence in the far zone when Nx » 1, » 1/ ^ 1 and 

Dy « 1. 

Without the knowledge of the function it would not be possible to fully 
predict the outcomes of a two-pinhole experiment. Results described here 
should be considered as an illustration of our general theory that may or 
may not, depending on the case under study, have practical influence. It 
may, however, be the subject of experimental verification. 



3.4.1 Discussion 

Although the Xf is independent of Nx and Dx, its actual shape is determined 
by the presence of a large horizontal emittance. Let us show this fact. If 
both Nx,y 1 and D^^y ^ 1 (filament beam limit), one would have had 

|G| = 01 02 j, so that \g\ = 1, strictly Note that in this case G 

could not be factorized. When instead Nx » 1 and » 1 the cross-spectral 
density can be factorized according to Eq. (56), but the integration in d(px, that 
follows from an integration over the horizontal electron beam distribution, 
is still included in the vertical cross-spectral density Gx (and, therefore, in 
Xp). This results in the outcome described above and can be traced back 
to the non-Gaussian nature of Wf. Note that if one adopted a Gaussian- 
Schell model, the cross-spectral density could have been split in the product 
of Gaussian intensity and Gaussian spectral degree of coherence and ^f, 
being a product of Gaussians, would have been separable. Then, \gy\ would 
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Fig. 16. Comparison between some zeros of Xv, at coordinates {dy,Ady (black 
circles), and the directivity diagram of undulator radiation in the vertical direction 
at very large horizontal electron beam divergence » 1 and negligible vertical 
divergence Dy <sc 1 (solid line). 

have been constant for Nx » 1 and D;^ » 1. As a result, the effect described 
here would not have been recognized. This fact constitutes a particular 
realization of our general remarks about Gaussian-Schell models at the end 
of Section 2.5. 

As a final note to the entire Section, we stress the fact that our theory of 
partial coherence in third generation light sources is valid under several non- 
restrictive assumptions. Alongside with previously discussed conditions 
» 1, Na, » 1, OtCO » N^u, and the assumption of perfect resonance (i.e. the 
limit C <?c 1), we assumed separability and particular shape of the electron 
beam phase space (see Eq. (52)). Moreover, for third generation light sources 
we assumed Cx ^ A (up to the VUV range). Together with (ix,y - l^w this 
implies Nx^ 1 and Dx^ 1, allowing for separability of the cross-spectral 
density in horizontal and vertical factors. Moreover, due to N^; » 1 and 
Dx » 1, we are dealing with quasi-homogenous Gaussian sources in the 
horizontal direction. This particular kind of sources will be treated in detail 
as an asymptote of our general theory in the next Section. In the vertical 
direction instead, we still have fully generic sources. We showed how the 
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vertical cross-spectral density can be expressed in terms of convolutions 
between two-dimensional universal functions and Gaussian functions. A 
particularly interesting case is that of quasi-homogeneous non-Gaussian 
sources that will also be treated as an asymptotic case in the following 
Section dedicated to quasi-homogeneous sources. 



4 Quasi-homogeneous asymptotes for undulator sources 

In Section 3 we developed a general theory of transverse coherence prop- 
erties of third generation light sources. In this Section we consider the class 
of quasi-homogeneous sources for undulator devices as an asymptotic limit 
for that theory. 

Quasi-homogeneous sources are defined by the fact that the cross-spectral 
density of the virtual source (positioned at z = 0) can be written as: 

G (O, r, Af) = 7 (O, r) g (Ar) . (80) 

The definition of quasi-homogeneity amounts to a factorization of the cross- 
spectral density as the product of the field intensity distribution and the 
spectral degree of coherence. A set of necessary and sufficient conditions 
for such factorization to be possible follows: (i) the radiation intensity at the 
virtual source varies very slowly with the position across the source on the 
scale of the field correlation length and (ii) the spectral degree of coherence 
depends on the positions across the source only through the difference Ar. 

Factorization of Eq. (55) as in Eq. (80), for third generation light sources, is 
equivalent to a particular choice of the region of parameters for the electron 
beam: Nj » 1, » 1 and either (or both) » 1 and Dy » In this 
case, the reader may verify that conditions (i) and (ii) are satisfied. 

In the horizontal direction, we have both Nx » 1 and » 1 for wave- 
lengths up to the VUV range, so that factorization of the cross-spectral 
density in horizontal and vertical contributions Gx and Gy always holds. 

Let us first consider Gy. Depending on the values of Ny and Dy we may 
have Gaussian quasi-homogeneous sources characterized by a Gaussian 

These conditions describe the totality of third generation quasi-homogeneous 
sources. In fact, while a purely mathematical analysis indicates that factorization 
of Eq. (55) is equivalent to more generic conditions (Nx » 1 and Ny » 1, or Dj » 1 
and Dy » 1), comparison with third generation source parameters (Nx » 1 and 
Dx » 1) reduces such conditions to the already mentioned ones. 
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transverse distribution of intensity {Ny » 1 and Dy » 1), as well as non- 
Gaussian quasi-homogeneous sources (Ny » 1 and Dy < 1 or Ny < 1 and 
Dy » 1). Gaussian quasi-homogenous sources are to be expected in the 
vertical direction in the hard X-ray limit, where diffraction effects play no 
role. On the contrary, diffraction effects must be accounted for when dealing 
with non-Gaussian quasi-homogeneous sources. 

Let us now consider Gx- In the horizontal direction both Nx ^ 1 and Dj; » 1. 
It follows that Gaussian quasi-homogeneous sources find a very natural 
application in the description of the cross-spectral density in the horizontal 
direction, from the hard X-ray to the VUV range. 

We will see that the VCZ theorem applies to all quasi-homogeneous cases. 
Actually, the concept of far-zone for quasi-homogeneous sources can be 
introduced as the region in the parameter space z, Nx,y, Dx,y such that the VCZ 
theorem holds. We will see that these condition coincides with condition (63) 
given before. 

4.2 Gaussian undulator sources 

When Nx » 1 and D;^ » 1 Eq. (57) applies. When Ny » 1 and Dy » 1, Eq. 
(58) reduces to 



that is equivalent to Eq. (57): thus, identical treatments hold separately in 
the horizontal and vertical directions. For this reasons, and for simplicity 
of notation, we drop all subscripts "x" or "y" in the present Section 4.1, 
and we substitute letters "x" and "y" in variables with the more generic 
"r". However, as stated before, the Gaussian quasi-homogeneous model 
primarily describes third generation light sources in the horizontal direction. 

Eq. (80) is satisfied. Moreover, 




(81) 




(82) 



and 



g{Ar) 



(Ar) 



(83) 



exp 



2 
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Propagation of Eq. (81) can be found taking the limit of Eq. (62) f or N » 1 
and D » 1, which yields an analytical expression for the evolution of the 
cross-spectral density based on Eq. (44) at C <^ 1 : 



G{z, r, Ar) = 



zVaTd 

xexp 



exp 

. ArAr 
'z(A + D) 



2(A + D)z2 



exp 



z J 



exp 



ADjArf ' 
2(A + D) 



(84) 



where 



(85) 



We have 



I{z,r) = 



zyfA + D 



exp 



2(A + D)z2 



(86) 



and 



g{z, Ar) = exp 



jAr 



exp 



. ArAr 
'z(A + D) 



exp 



AD(Ar)2 



2(A + D) 



(87) 



Note that, due to the phase factors in Eq. (84), only the virtual source at z = 
constitutes a quasi-homogeneous source. 

Geometrical interpretation of A is the dimensionless square of the apparent 
angular size of the electron bunch at the observer point position, calculated 
as if the beam was positioned at z = 0. 

The far-zone for quasi-homogeneous Gaussian sources is given by condition 
A «: D, which can be retrieved by condition (63) or directly by Eq. (84). In 
this case, simplification of Eq. (84) or use of Eq. (65) in the limit f or N » 1 
and D » 1 yields the far-zone cross-spectral density: 



G (z, 0, Ad) = I exp [izOAO] exp 



N(A0)2" 






2 


exp 


2D 



(88) 



so that 
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and 



N{Ad)' 
2 



,21 



g (z. Ad) = exp [izdAd] exp 



(90) 



Similarly as before we suppressed subscripts "x" or "y" in the symbol 6. 

Analysis of Eq. (81) and Eq. (88) allows to conclude that 

(^2) the spectral degree of coherence of the field at the source plane g{Ar) 
and the angular distribution of the radiant intensity 1(6) are a Fourier pair. 

(&2) the spectral degree of coherence of the far field g{Ad) and the source- 
intensity distribution I{f) are, apart for a simple geometrical phase factor, a 
Fourier pair. 

The statement (^2) is a version of the VCZ theorem valid for quasi-homogeneous 
sources. Statement (fl2) instead, regards the symmetry between space and 
angle domains, and can be seen as an inverse VCZ theorem. 

This discussion underlines the link between the VCZ theorem and the 
Wiener-Khinchin theorem. Exactly as the space domain has a reciprocal 
description in terms of transverse (two-dimensional) wave vectors, the time 
domain has a reciprocal description in terms of frequency. The reader will 
recognize the analogy between statements (^2) and (^2)/ with statements (fli) 
and (bi) discussed in Section 2.3. In particular, the VCZ is analogous to the 
inverse Wiener-Khincin theorem. Similarly, separability of G in Eq. (80) in 
the product of spectral degree of coherence and intensity is analogous to 
separability of r^,,, in Eq. (11) in the product of spectral correlation function 
and spectral density distribution of the source. 

Let us calculate the transverse coherence length as a function of the 
observation distance z. We introduce the coherence length following the 
definition by Mandel [14]. The coherence length, naturally normalized to 
the diffraction length ylL^c/o) is defined as 



00 




(91) 



—00 
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Fig. 1 7. Coherence length tc as a function of z and asymptotic behaviors for z — > 1/2 
and z » 1. Here N = 10^ and D = 10. 

Performing the integration in Eq. (91) with the help of Eq. (87) yields: 
4(z) = V^(x + ^) • (92) 



The coherence length in Eq. (92) exhibits linear dependence on z, that is 
— ^ yfnjN z while for z — > 1/2 that is at the end of the undulator, it 
converges to a constant tc — ^ [tt/ (4N) + tt/D]^/^. Eq. (92) and its asymptotes 
are presented in Fig. 17 for the case N = 10^, D = 10. At the exit of the 
undulator, ~ 1/ Vd, because N » D. On the other hand, horizontal 
dimension of the light spot is simply proportional to Vn. This means that 
the horizontal dimension of the light spot is determined by the electron beam 
size, as is intuitive, while the beam angular distribution is printed in the fine 
structures of the intensity function that are of the dimension of the coherence 
length. In the limit for A D the situation is reversed. The radiation 
field at the source can be presented as a superposition of plane waves, 
all at the same frequency a», but with different propagation angles with 
respect to the z-direction. Since the radiation at the exit of the undulator is 
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Fig. 18. Physical interpretation of the generalized VCZ theorem. If the radiation 
beyond the source plane is partially coherent, a spiky angular distribution of inten- 
sity is expected. The nature of these spikes is easily described in Fourier transform 
notations. We can expect that typical width of the angular distribution of intensity 
should be of order (cuA/c)"^, where A is the typical linear dimension of spatially 
random intensity fluctuations. If the source has transverse size d, the angular dis- 
tribution of intensity should contain spikes with typical width of about {codjcY^, a 
consequence of the reciprocal width relations of Fourier transform pairs. 

partially coherent, a spiky angular distribution of intensity is to be expected. 
The nature of the spikes is easily described in terms of Fourier transform 
theory. From Fourier transform theorem we can expect an angular spectrum 
with Gaussian envelope and rms width Vd. Also, the angular distribution 
of intensity should contain spikes with characteristic width 1/ Vn, as a 
consequence of the reciprocal width relations of Fourier transform pairs 
(see Fig. 18). This is realized in mathematics by the expression for the cross- 
spectral density, Eq. (84) and by the equation for the coherence length, Eq. 
(92). 

It is also important to remark that the asymptotic behavior for A «: 1 of g, 
that is Eq. (90) and 4, that is 




are direct application of VCZ theorem. In fact, the last exponential factor on 
the right hand side of Eq. (84) is simply linked with the Fourier transform 
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of // [Ij. We derived Eq. (84) for N » 1 and D » 1, with A «: D: ir\ non- 
normalized units these conditions mean that the VCZ theorem is applicable 
when the electron beam divergence is much larger than the diffraction angle, 
i.e. o'^ » A/(2nL,,,), the electron beam dimensions are much larger than the 
diffraction siz^i^ i.e. » AL,^lln, and {o'zf » cr^. 

In list] (paragraph 5.6.4) a rule of thumb is given for the applicability region 
of the generalization of the VCZ theorem to quasi-homogeneous sources. 
The rule of thumb requires z > 2rfA/A where d is "the maximum linear 
dimension of the source", that is the diameter of a source with uniform 
intensity and A "represents the maximum linear dimension of a coherence 
area of the source". In our case, since a is the rms source dimension, d — 2o. 
Moreover, from Eq. (92) we have A = - A/(2 ^Jno'). The rule of thumb 
then requires z > 2cr/( ^fno'): in dimensionless form this reads z > -\/NJD. 
This is parametrically in agreement with our limiting condition A « D, even 
though these two conditions are different when it come to actual estimations: 
our condition A « D is, in fact, only an asymptotic one. To see how well it 
works in reality we might consider the plot in Fig. 17. There N = 10^ and D = 
10. Following [3] we may conclude that a good condition for the applicability 
of the VCZ theorem should be z > 10. However as it is seen from the figure, 
the linear asymptotic behavior is not yet a good approximation at z ^ 10. 
This may be ascribed to the fact that the derivation in |j^ is not generally 
valid, but has been carried out for sources which drop to zero very rapidly 
outside the maximum linear dimension d and whose correlation function 
also drops rapidly to zero very rapidly outside maximum linear dimension 
A. However, at least parametrically, the applicability of the VCZ theorem 
in the asymptotic limit A «: D can be also expected from the condition 
z > 2rfA/A in [3]. 

We conclude that the far field limit A «: D corresponds with the applicability 
region of the VCZ theorem. When this is the case, the VCZ theorem applies 
and the modulus of the spectral degree of coherence in the far field, forms a 
Fourier pair with the intensity distribution of the virtual source. In particular 
one concludes that the rms width of the virtual source is Vn. In our study 
case for N » 1 and D » 1, such a relation between the rms width of the 
spectral degree of coherence in the far field and the rms dimension of the 
virtual source is also a relation between the rms width of the cross-spectral 
density function in the far field and the rms dimension of the electron beam 
at the plane of minimal beta function in the center of the undulator. In 
dimensional units one can write the value Oc of the rms width of the spectral 

We do not agree with statement in |@]: "the electron-beam divergence must be 
much smaller than the photon divergence" for the VCZ theorem to apply- This 
would imply o' «; -^jA/ilnL^) (reference [7^, page 571, Eq. (57)). 
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degree of coherence g{Ar) in the far field as 



Here o is, as usual, the rms dimension of the electron beam. These few last 
remarks help to clarify what is the size of the source in the VCZ theorem. 



that is far from being a trivial question. In several papers |19l, |2O0 the rms 
electron beam size is recovered from the measurement of the transverse 
coherence length and subsequent application of the VCZ theorem, under the 
assumptions that the VCZ theorem can indeed be applied. In this regard, in 
| }l9l ] Section V, one may find a statement according to which the rms electron 
beam size "is only the average value along the undulator" because "the 
beta function has a large variation along the undulator". Another example 



dealing with the same issue is given in reference [2011 . This paper (as well 
as reference ||l9ll) reports experimental results. However, authors of | |2otl 
observe a disagreement between the electron beam rms size reconstructed 
from the VCZ theorem and beam diagnostics result of about a factor 2. They 
ascribe this variation to the variation of the electron beam size along the 
undulator. In footnote 25 of reference [20], one may read: "The precise shape 
and width of the x-ray intensity distribution in the source plane are directly 
connected to the properties of the electron beam. It would not be surprising if 
the limited depth of focus of the parabolically shaped electron beta function 
in the undulator translates into a virtually enlarged x-ray source size.". At 
first glance it looks like if the SR source had a finite longitudinal dimension, 
and the virtual source size depended on variations of the beta function 
along the undulator. However, as we have seen before, the concept of virtual 
source involves a single transverse plane, and in the most general case any 
variation of the beta function does not affect the virtual source size. In our 
case of quasi-homogeneous Gaussian source, the virtual source is located 
where the beta function is minimal, and its size coincides with the transverse 
size of the electron beam at that location. 

Finally, it should be remarked that the VCZ theorem can only be applied to 
quasi-homogeneous sources. Consider for example 21], where a char- 



acterization of the vertical emittance in Spring-8 is reported, based on the 
measurement of the X-Ray beam coherence length in the far zone. The ex- 
periment was performed at the beamline BL29XU. Based on the assumption 
of validity of the VCZ theorem, it was found that the rms electron beam size 
at the undulator center (corresponding to the minimal value of the beta 
function) was Sy - 4.5 /im, and that the coupling factor between horizontal 
and vertical emittance was down to the value C - 0.12%, which corresponds 
to an extremely small vertical emittance Cy = 3.6 pm-rad. A resolution limit 
of this method was also discussed, based on numerical calculations of the 
radiation size from a single electron performed at the exit of the undulator. 



53 



Sp — 1.6 [im at Ep = 14.41 keV for the 4.5 m long undulator used in the 
experiment. The resolution limit of the measurement of Sy was estimated to 
be about 1 /.tm. Considering propagation laws for Gaussian beams it seems 
reasonable that the radiation spot size at the virtual source, located in the 
center of the undulator, be smaller than Sp. However, undulator radiation 
from a single electron cannot be considered a Gaussian beam. In particular, 
under the resonance approximation the field at the exit of the undulator 
exhibits a singularity and is not suitable for evaluation of the radiation spot 
size. 

Use of Sp led to an underestimation of the virtual source size and, thus, 
an overestimation of the resolution. Let us show this fact. Based on Eq. 
(30), we can determine the correct virtual source size of undulator radiation 
from a single electron. Let us consider the case when the single electron 
is emitting photons at the fundamental harmonic with energy E = 14.41 
keV. The angular frequency of light oscillations is given, in this case, by 
CO = 2.2 ■ 10^^ Hz. For an undulator length L;„ = 4.5 m, the normalization 
factor for the transverse size, (LiyC/o;)^^^, is about 8 /.im. From Fig. 8 obtain 
the dimensionless Half Width Half Maximum (HWHM) radiation size from 
a single electron (i.e. the HWHM width of the intensity distribution at 
the virtual source, located at the center of the undulator). This HWHM 
dimensionless value is about 0.7. It follows that the HWHM value of the 
radiation spot size from a single electron is about 0.7 • {cLu,/ojY^^ - 6 /.im. 
Therefore, the resolution in [|19^21t]/ is estimated to be better than the correct 
value. 



Note that the HWHM radiation spot size from a single electron is larger 
than the rms electron beam size Sy — 4.5 y.m. found by means of coherence 
measurements. One concludes that the uncertainty due to finite resolution is 
larger than the measured electron beam size. This suggests that the analysis 
of experimental results |19, 2l[] might include a logical flaw. Authors of that 
reference assume the validity of the VCZ theorem in the vertical direction. 
If one assumes their result of a vertical emittance Cy — 0.3A/(27i), it follows a 
posteriori that the VCZ theorem could not have been applied in first instance 
(in this experiment the value of the beta function was jS ^ Lj„). We suggest 
that analysis of experimental results in |19l,l2ltl should be based, instead, on 
the study of transverse coherence for non-homogeneous undulator sources 
in free space made in the previous Section 3. 



4.2 Non-Gaussian undulator sources 



Let us now turn to the analysis of non-Gaussian quasi-homogeneous sources 
on the basis of Eq. (58) and Eq. (65). We still assume that Nx:$>l and Dy^:$>l. 
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Moreover, as before, we assume that the minimal beta function is located at 
the undulator center. 



4.2.1 Case of a large electron-beam size Ny » 1 and Dy < 1. 
Eq. (58) simplifies to 
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(95) 



where we defined the universal function ysioi) as 
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(96) 



A plot of ysicx) is given in Fig. 19. Main features of ys are a strong non- 
Gaussian shape, and the fact that also negative values are assumed. Eq. (80) 
is satisfied. Moreover, 
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(97) 



and 



|y(0,Ay) = exp 



(Ay)^D, 



7s(Ay) , 



(98) 



In the far-zone, Eq. (65) simplifies to 



Gy(Z,0y,A0y) = 



exp 



exp 
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a 



Fig. 19. Plot of the universal function ys{oi), used to calculate the cross-spectral 
density of a quasi-homogeneous source when » 1, D;,: » 1/ ^'^y » 1 and Dy < 1. 
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(99) 



where Jp has already been defined in Eq. (13). Moreover, 



Inz- 



00 



(0y + ^y)^ 



2D. 



^F((^)y) , 



(100) 



and 



gy{Z, t^Qy) = [zZ0yA0y] CXp 



NyA02 



(101) 



Note that, due to geometrical phase factor in Eq. (99), only the virtual source 
at z = constitutes a quasi-homogeneous source. 
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Also, ys and If basically form a Fourier pair: 



00 



(102) 



meaning that the inverse VCZ theorem is satisfied for Nx,y <c z^. 

Finally, it is possible to calculate ys analytically. To this purpose, it is suffi- 
cient to note that the Fourier transform 
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-00 —00 



can be evaluated with the help of the Bessel-Fourier formula as 
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where A' = + rf, cj)^ = (pl + (p^, Si(-) is the sine integral function and Ci(-) 
is the cosine integral function. Thus, letting = and r/ = Ay one has 



7s(Ay) = - 
n 



n + (Ay)^Ci 
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2Si 



./(Ay)' 



(105) 



4.2.2 Case of a large electron-beam divergence Dy » 1 and Ny < 1. 

A similar analysis can be given in the case of a large electron-beam diver- 
gence. In this case Eq. (58) simplifies to give 
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(106) 



where Is has already been defined in Eq. (70). Eq. (80) is satisfied. Moreover, 
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^y(0,Ay) = exp 



(108) 



In the far zone, Eq. (65) simplifies to 

G(z, 0y, A0y) = - Jy^-exp [iz0yA0y] exp 
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A plot of yf (a) is given in Fig. 20. As for ys, main features of yp are a strong 
non-Gaussian shape, and the fact that also negative values are assumed. Due 
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Fig. 20. Plot of the universal function yp, used to calculate the cross-spectral density 
in the far zone when Nx » 1, » 1, Ny < 1 and Dy » 1. 



to the phase factors in Eq. (109), only the virtual source at z = constitutes 
a quasi-homogeneous source. 



Also, Is and yp form a Fourier pair: 



I. 



oo 

;(y) = j d(py exp [-iy(py] yF{(py) , 



(113) 



meaning that the VCZ theorem is satisfied for z^Dx,y » 1. Since Dx,y » 1 this 
means that, in this case, the far zone begins at the very exit of the undulator, 
atz ~ 1. 



Finally, starting with the representation of 1$ in Eq. (70), we can write (see 
10] Appendix A.3.) 
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oo oo 

Is(y) = -Ks exp [-iy^,] j ia ah(a<p,)%ia) 



(114) 



Comparison with Eq. (113) yields the following alternative representation 
of in terms of a one-dimensional integral involving special functions: 



oo 

yA(\>y) = n J da a/o(a(/)y)Wo(a) 



(115) 



4. 3 Accuracy of the quasi-homogeneous approximation 



As we have discussed before, the quasi-homogeneous approximation can 
be applied when either or both Ny » 1 or Dy » 1. In this Section we will 
see that the accuracy of the quasi-homogeneous approximation scales as 
( ^Jmax[Ny, 1] max[Dy, 1] )"^. As we have previously seen, in the particular 
case when both Ny » 1 and Dy » 1, a Gaussian model may be used. We 
will see that the accuracy of such model is worse than that of the quasi- 
homogeneous model, and scales as max(l/ JD^,1/ jNy). 



We begin demonstrating that, when the source is quasi-homogeneous, we 
may take the approximation S{Ay,^) - 'ys{Ay)Is{^) in Eq. (60) with an 
accuracy scaling as ( ymax[Ny, 1] max[Dy, 1] )"^. First, let us introduce a 
normalized version of the one-dimensional inverse Fourier transform of the 
function 'F, that is 
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where the normalization factor is chosen in such as way that !F'(0,0) = 1. 
The cross-spectral density in Eq. (60) can therefore be written as 
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having used the convolution theorem. Under the quasi-homogeneous as- 
sumption, we can approximate !F(Ai/, u) - !F(Ai/, 0)9^(0, u). To show this, let 
us represent 'F{u,v) using a Taylor expansion around the point (0,0). One 
obtains 



1 

k=l 

+0{uv) , 



+ v 



u=0 



(118) 



where the normalization relation ;F(0, 0) = 1 has been taken advantage of. 
Similarly, one may consider the following representation of the product 
!F(m, 0) 1^(0, V) also obtained by means of a Taylor expansion: 
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(119) 



having used !F(0,0) = 1. Comparison of the last equality in (119) with 
the right hand side of Eq. (118) shows that 'f'{u,v) ^ 'F{u,0)T'{0,v) up 
to corrections of order uv ~ 1/ -^max[Ny, 1] max[Dy, 1], that is the quasi- 
homogeneous accuracy. Using this approximation in Eq. (117) yields 
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Finally, recalling the definitions of ys and 1$ we can write Eq. (120) as 
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Eq. (121) is valid in any quasi-homogeneous case. 

Note that Eq. (121) accounts for diffraction effects through the universal 
functions ys and Is- This may be traced back to the use of the inhomoge- 
neous wave equation to calculate the cross-spectral density for the virtual 

source, from which Eq. (121) follows. Deriving Eq. (121), we assume a large 
number of modes, and this justifies the use of phase space representation as 
an alternative characterization of the source, in place of the cross-spectral 
density (i.e. Eq. (121) itself). 

Setting Ay = 0, Eq. (121) gives the exact intensity distribution at the virtual 
source, i.e. Eq. (68). The spectral degree of coherence on the virtual source 
is then recovered using the definition of quasi-homogeneous source G = 
I{y)g{Ay). Since the source is quasi-homogeneous, the Fourier transform of 
the spectral degree of coherence g{Ay) yields the intensity in the far zone. 
Remembering that ys and Ip form a Fourier pair, we conclude that, starting 
from Eq. (121) it is possible to reproduce the exact result for the intensity 
in the far zone, Eq. (71). Quite remarkably, Eq. (121), which is derived 
under the quasi-homogeneous approximation and is related to an accuracy 
1/ ymax[Ny, 1] max [Dy, 1], yields back two results, Eq. (68) and Eq. (71) 
which are valid regardless the fact that the source is quasi-homogeneous or 
not. 

Let us now consider the quasi-homogeneous case when both Ny » 1 and 
Dy » 1. Worsening the accuracy in the calculation of the cross-spectral 
density of the source, we may reduce Eq. (121) to 



(5y(0,y,Ay)= -^^exp 

that is Eq. (81). Note that neglecting the product with the ys function can be 
done with an accuracy 1 / ^fD^, while extraction of the exponential function 
in y from the convolution product with the Is function can be done with an 
accuracy 1/ ^/N^■ In our study case when Dy » 1 and Ny » 1, the overall 
accuracy of Eq. (122) can be estimated as max(l/ yjD^,l/ ^^), that is the 
accuracy of the Gaussian approximation. Such accuracy is much worse than 
that of the quasi-homogeneous assumption in Eq. (121), that is 1/ -^N^Dy. 

When Ny » 1 and Dy - 1 the accuracy of the quasi-homogeneous ap- 
proximation becomes 1 / ^JNymax[l,Dy]. When Ny ^ 1 and Dy » 1 it be- 
comes, instead, 1/ -^max[l, Ny]Dy. In these cases, the accuracy of the quasi- 
homogeneous approximation is comparable to the accuracy of the Gaussian 
approximation. To be specific, when Ny » 1 and Dy =^ 1 the accuracy of the 
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quasi-homogeneous approximation is of order 1/ i^Ny and Eq. (121) can be 
substituted with 















exp 


2 



ys(Ay) 



(123) 



without loss of accuracy, because the relative accuracy of the convolution is 
of order 1 / -^JN^ as the accuracy of the quasi-homogenous approximation. 
Eq. (123) is just Eq. (95). A similar reasoning can be done when Dy » 1 and 
Ny ^ 1. In this case the accuracy of the quasi-homogeneous approximation 
is of order 1/ JD^, and Eq. (121) can be substituted with 



G(0,y,Ay) = 
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Isicp) (124) 



without loss of accuracy. Eq. (124) is just Eq. (106). 

4.4 Quasi-homogeneous sources in terms of phase space 

The cross-spectral density at the virtual source can be written as in Eq. (80) 
that we rewrite here for convenience in terms of coordinates fx,y and Arxy. 



Goirx, ry, Arx, Ary) = I (r^, fy) giAr^, Ary) . 



(125) 



For notational simplicity we use notation Go to indicate G at z = 0. The 
Fourier transform of Eq. (125) with respect to all variables will be indicated 
with 



oo oo 



go{9x,9y,A9x,A9y) = j dAr', j dAr'y j df, j df'y Go{f'„f'y, Ar'„ Ar'y 



) 



-oo — oo 



X exp[2z(0,Ar; + ByAr'y)] exp[i(A0,f; + A0yf;)] . 

(126) 

The two quantities l{rx, fy) = Go{fx, fy, 0, 0) and t{9x, 9y) = Qo{9x> 9y, 0, 0) are 
always positive because, by definition of Go, they are ensemble averages of 
quantities under square modulus. 



Let us now introduce the Fourier transform of Eq. (80) with respect to Afx^y: 
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— oo — oo 

xexp[i{dAK + eyAr'y)]. (127) 

Accounting for Eq. (125), i.e. in the particular case of a (virtual) quasi- 
homogeneous source, Eq. (127) can be written as 



Oo(f., fy, 0„ By) = I (f,, fy) f (0„ By) , (128) 

having recognized that t{Bx, By) = Qo{Bx, By, 0, 0) is the Fourier transform of 
the spectral degree of coherence g. The distribution Oo, being the product 
of two positive quantities, never assumes negative values. Therefore it may 
always be interpreted as a phase space distributiorl^. This analysis shows 
that quasi-homogeneous sources can always be characterized in terms of 
Geometrical Optics. It also shows that, in this particular case, the coordinates 
in the phase space, fx,y and Bx,y, are separable. 

Eq. (127) is the definition of a Wigner distribution. In the case of quasi- 
homogenous sources, as we have just seen, the Wigner distribution is never 
negative and, therefore, can always be interpreted as a phase space distri- 
bution. In the case of non quasi-homogeneous sources one may still define 
a Wigner distribution using Eq. (127). However the Wigner function itself is 
not always a positive function. As a consequence it cannot always be inter- 
preted as a phase space distribution. On the one hand, quasi-homogeneity 
is a sufficient condition for the Geometrical Optics approach to be possibly 
used in the representation of the source. On the other hand though, nec- 
essary and sufficient conditions for Oo to be a positive function are more 
difficult to find. 

Note that in the case of quasi-homogeneous non-Gaussian sources one 
should account to diffraction effects when calculating source properties. 
However, using a Wigner function approach, one can still use a phase-space 
representation. In this case diffraction effects have the effect of complicat- 
ing the structure of the phase-space describing the source. Also note that 
Gaussian-Schell model cannot be applied in all generality even to quasi- 
homogeneous sources, as it does not properly describe the case of non- 



Physically, in the quasi-homogeneous case, f can be identified with the radiant 
intensity of the virtual source. This follows from a statement similar to the van 
Cittert-Zernike theorem for quasi-homogeneous sources (see ^). Note that the 
intensity and the Fourier transform of the spectral degree of coherence are obtained 
back from the phase space distribution, Eq. (128), by integration over coordinates 
6x,y and fx,y respectively. 
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Gaussian sources, which is most natural for third generation facilities in the 
vertical direction. 

As we have seen before, third generation light sources are characterized, up 
to the VUV range, hyNx^l and » 1, which implies factorization of the 
cross-spectral density. Thus, considering the vertical direction separately, 
equivalent condition for quasi-homogeneity requires that either (or both) 
Ny » 1 and Dy » 1. 

An intuitive picture in the real space is given by a (virtual) quasi-homogeneous 
source with characteristic (normalized) square sizes of order max[Ny, 1] and 
Nx, and characteristic (normalized) correlation length square sizes of order 
min[l/Dy, 1] and l/D^, in the vertical and horizontal directions. Since hori- 
zontal and the vertical directions can be treated separately, we have a large 
number of independently radiating sources given by the product 

My = max[Ny, 1] max[Dy, 1] (129) 

in the vertical direction, and 

Mx = NxDx (130) 

in the horizontal direction. The number Mx^y is, in other words, an estima- 
tion of the number of coherent modes in the horizontal and in the vertical 
directiorl^ The number M^y is the accuracy of Geometrical Optics re- 
sults compared with Statistical Optics results or, better, the accuracy of the 
quasi-homogeneous assumption. It should be noted that, as Mx,y approaches 
unity, the accuracy of the quasi-homogeneous assumption becomes worse 
and worse and Mx,y cannot be taken anymore as a meaningful estimation 
of the number of modes: it should be replaced by a more accurate concept 
based on Statistical Optics. To complete the previous statement we should 
add that Mx,y completely loses the meaning of "number of modes" when 
Geometrical Optics cannot be applied. For instance when both Ny and Dy 
are of order unity (or smaller), one can state that the Geometrical Optics 



This is in agreement with an intuitive picture where the photon-beam phase 
space reproduces the electron-beam phase space up to the limit imposed by the 
intrinsic diffraction of undulator radiation. Imagine to start from a situation with 
Nx,y » 1 and Dx,y » 1 and to "squeeze" the electron-beam phase space in the 
vertical direction by diminishing Ny and Dy. On the one hand the characteristic 
sizes of the phase space of the electron beam are always of order Ny and Dy. On 
the other hand the characteristic sizes of the phase space of the photon beam are 
of order max[Ny, 1] and max[Dy, 1]: diffraction effects limit the "squeezing" of the 
phase space of the photon beam. 
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approach fails in the vertical direction because the phase space area is get- 
ting near to the uncertainty limit. In this case it is not possible to ascribe 
the meaning of "number of modes" to the number My simply because the 
Geometrical Optics approach in the vertical direction fails. However, when 
Ny and Dy are of order unity (or smaller), but both 1 and » 1, the 
cross-spectral density admits factorization in the horizontal and in the verti- 
cal direction and the source in the horizontal direction can be still described, 
independently, with the help of Geometrical Optics. 

We should remark that Statistical Optics is the only mean to deal, in general, 
with the stochastic nature of SR. Only in particular cases SR can be treated 
in terms of Geometrical Optics, where planning of experiments can take 
advantage of ray-tracing code techniques. One of these cases is constituted 
by second generation light sources, because Nx,y » 1 and Dx,y » 1. 

Description in terms of quasi-homogeneous sources is also important for 
third generation facilities in the case of bending magnet beamlines (surpris- 
ingly, in both horizontal and vertical plane), horizontal plane in undulator 
beamlines, but it cannot be applied with accuracy to future sources (e.g. 
ERL-based sources). 

Quasi-homogenous sources can be described in terms of geometrical optics, 
the outcome being equivalent to description in terms of statistical optics. In 
order to decide whether Geometrical Optics or Wave Optics is applicable, 
in all generality one should separately compare the photon beam size and di- 
vergence with the radiation diffraction size and diffraction angle, which are 
quantities pertaining the single electron radiation. Let us fix a given direc- 
tion X or y. The square of the diffraction angle is defined by (c'^)^ ~ A/{2nLf), 
Lf being the formation length of the radiation at wavelength A. The diffrac- 
tion size of the source is given by oa ~ o'^Lf. In calculating the photon 
beam size and divergence one should always include diffraction effects. 
As a result, if and (cr')^ indicate the square of the electron beam size 
and divergence, the corresponding square of the photon beam size and 
divergence will be respectively of order max[cr^,a^] and max[(cr')^, (cr^)^]. 
These quantities can be rewritten in terms of the electron beam emittance 
as max[e:j3, o^] and max[e/ jS, (cr^)^], jS being the beta function value at the vir- 
tual source position for the radiator (undulator, bending magnet, or other). 
Dividing these two quantities respectively by and {o'^)^ give natural val- 
ues, normalized to unity, for the photon beam size max[2ne(i/{LfA), 1] and 
divergence max[2HeL^/(j3A), 1]. When the product between these two quan- 
tities is much larger than unity one can use a Geometrical Optics approach. 
In this case, this product represents the normalized photon beam emittance. 
When j8 ~ L^^, as in many undulator cases, one may compare, for rough 
estimations, the electron beam emittance and the radiation wavelength as 
we have done before. However, in the case of a bending magnet one may 



66 



typically have jS of order 10 m and Lf - (p^A)^^^ (p being the bending 
radius) of order lO"'^ -r 10"^ m. The ratio (}/Lf » 1 now constitutes an ex- 
tra large parameter of the problem. In this case, even if the electron beam 
emittance is two order of magnitude smaller than the wavelength, due to 
diffraction effects one can still apply a Geometrical Optics approach, because 
max[27iejS/(L|A), l]-max[2neLf/{^A), 1] » 1, i.e. the photon beam emittance 
is much larger than the wavelength. As a result, dimensional analysis sug- 
gests that bending magnet radiation may be treated exhaustively in the 
framework of Geometrical Optics even for third generation light sources. 



5 Conclusions 

This work presents a theory of transverse coherence properties from third 
generation light sources, valid while radiation evolves in free-space. Besides 
being important for experiments involving coherence that make no use of 
optical elements, it constitutes the first step towards the solution of the 
image formation problem for undulator sources (see [22]) that will be a 
natural follow-up to the present article. 

We considered Synchrotron Radiation (SR) as a random statistical process 
to be described using the language of statistical optics. Statistical optics 
developed in connection with Gaussian, stationary processes characterized 
by homogeneous sources. However, for SR, there is no a priori reason to 
hold these assumptions satisfied. 

We showed that SR is a Gaussian random process. As a result, statistical 
properties of SR are described satisfactory by second-order field correlation 
functions. We focused, in particular, on undulator sources. It should be 
noted here that wiggler and bending magnets are still being used at third 
generation facilities. However, as has been remarked in the previous Section 
4.4, these devices are characterized by a much shorter formation length, 
which allows one to apply a Geometrical Optics approach to describe them. 
Thus, in this case, the formalism developed for second generation facilities 
can be satisfactory taken advantage of. In other words, use of wigglers 
and bending magnets is mainly related with applications requiring higher 
photon flux (but not high coherent flux), compared with analogous devices 
installed in second generation facilities. Our choice of considering undulator 
sources is justified by the fact that we are interested in transverse coherence 
properties of radiation. 

With this in mind, a frequency-domain analysis was used to describe un- 
dulator sources from a mathematical viewpoint. As a consequence of the 
frequency domain analysis we could study the spatial correlation for a given 
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frequency content using the cross-spectral density of the system. This can 
be used to extract information even if the process is non-stationary, and 
independently of the spectral correlation function. 

We gave a general expression for the cross-spectral density dependent on 
six dimensionless parameters. Subsequently we assumed small normalized 
detuning from resonance, thus obtaining a simplified expression of practical 
relevance. 



We simplified our expressions further in the case of third generation light 
sources, based on a large horizontal emittance (compared with the wave- 
length). In this case, the cross-spectral density can be factored in the product 
of a horizontal and a vertical factor. 

Attention was subsequently drawn on the vertical cross-spectral density Gy, 
without loss of generality. We expressed Gy in terms of one-dimensional con- 
volutions between universal functions and Gaussian functions, and studied 
different asymptotic cases of interest. 

In the case of a vertical emittance much smaller than the radiation wave- 
length we derived the counter-intuitive result that radiation is not fully 
coherent in the vertical direction. This effect can be interpreted as an influ- 
ence of a large horizontal emittance on the vertical plane, and is related with 
the particular non-Gaussian nature of the single particle field. 

Subsequently, we studied quasi-homogeneous cases of interest, discussing 
the applicability of the VCZ theorem for Gaussian and non-Gaussian quasi- 
homogeneous sources. Finally, we discussed the accuracy of the quasi- 
homogeneous model. 

It is interesting to spend a few words about relation of our theory with nu- 
merical techniques. Computer codes have been written (see e.g. |23 , 23]) that 



deal with beamline design, based on wave-optics techniques. Codes also be- 
gin to be used to treat the case of partially coherent radiation: one of their 
final goals is to solve the image formation problem starting from first prin- 
ciples. Results may be obtained using numerical techniques alone, starting 
from the Lienard-Wiechert expressions for the electromagnetic field and ap- 
plying the definition of the field correlation function without any analytical 
work. Calculation of the intensity at a single point on the image plane can 
be approached by propagating wavefronts from different macro-particles 
through the entire optical system, calculating intensities and summing them 
up. This method is well-suited for parallel processing, and relatively easy to 
implement. Yet, a first-principle calculation of the field correlation function 
between two generic points at the image plane involves very complicated 
and time-expensive numerical evaluations. To be specific, one needs to per- 
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form two integrations along the undulator device and four integrations 
over the electron-beam phase space distribution to solve the problem in 
free space. Then, even in the simple case when the optical beamline can 
be modelled as a single focusing lens, other four integrations are needed 
to characterize coherence properties on the image plane, for a total of ten 
integrations. The development of a universal code for any experimental 
setup is then likely to be problematic. A more conservative approach may 
suggest the use of computer codes based on some analytical transformation 
of first principle equations suited for specific experimental setups. From 
this viewpoint our most general expressions (or alternatively, as it is being 
done, expressions for the Wigner distribution function) may be used as reli- 
able basis for the development of numerical methods. Our analytical theory 
allows treatment and physical understanding of many asymptotes of the 
parameter space and their applicability region with the help of a consistent 
use of dimensional analysis. This physical understanding, together with 
the possibility of using our asymptotic results as benchmarks for numerical 
methods, will be of help to code writers. 

In closing, as has been remarked in [6]: "[...] it is very desirable to have 
a way to model the performance of undulator beamlines with significant 
partial coherent effects, and such modelling would, naturally, start with 
the source. The calculation would involve the knowledge of the partial 
coherence properties of the source itself and of how to propagate partially 
coherent fields through space and through the optical components used in 
the beamline. [...] it is important to recognize that, although most of these 
calculations are, in principle, straightforward applications of conventional 
coherence theory (Born and Wolf, 1980; Goodman, 1985), there is not much 
current interest in the visible optics community". These statements were 
formulated more than ten years ago, when operation of third generation 
light sources started. While it was immediately recognized that usual SR 
theory was not adequate to describe them, no theoretical progress was 
made in that direction. Our paper finally answers the call in [6]. 
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